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We present a formalism and explicit results for two-loop flavor singlet evolution kernels of 
skewed parton distributions in the minimal subtraction scheme. This approach avoids explicit 
multiloop calculations in QCD and is based on the known pattern of conformal symmetry breaking 
in this scheme as well as constraints arising from the graded algebra of the TV = 1 super Yang- 
Mills theory. The conformal symmetry breaking part of the kernels is deduced from commutator 
relations between scale and special conformal anomalies while the symmetric piece is recovered 
from the next-to-leading order splitting functions and M = 1 supersymmetry relations. 
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1 Introduction 



A main source of information about the structure of matter is gained from reactions involving 
hadrons, with the theoretical description being most reliable for processes with one or more hard 
scales as for example the photon momentum transfer square, Q 2 , in ep reactions. The lowest 
approximation for a process, when only the Fock components with the minimal number of con- 
stituents in the hadron wave function at tree level are accounted for, corresponds to a simple 
picture of Feynman's parton model, which neglects QCD dynamics. Once this is taken into ac- 
count, it changes (sometimes drastically) the results. The specific predictions of QCD are given 
by scaling violations phenomena which arise due to two sources which play a complementary 
role depending on the magnitude of the hard momentum transfer. For moderately large Q 2 it is 
enough to include higher order perturbative corrections which scale as logQ 2 . On the other hand 
approaching the low Q 2 ~ few GeV 2 domain one has to take into account contributions relatively 
suppressed by powers of the momentum transfer 1/Q T ~ 2 , with r being the twist of the operators 
contributing to the amplitude of the process. Obviously, in this region this latter Q 2 behavior 
overwhelms the weaker log Q 2 due to perturbative evolution. However, even at rather small Q 2 the 
perturbative dependence on Q 2 cannot be discarded: being small at large Q 2 it becomes prominent 
at small Q 2 and modifies in a significant way the shape of the leading parton configuration of the 
hadronic wave function. Therefore, the knowledge of perturbative evolution is indispensable for 
the construction of the leading twist component from the genuinely non-perturbative distribution 
at low scales. Although this evolution can be weaker than the real non-perturbative evolution - 
which has to be used at very low Q 2 ~ "ihadron — ^ gives the right direction of change and can be 
seen as a builder of the rough features of hadronic distributions or amplitudes. The same applies 
for higher twists, but these contributions are usually discarded. 

In the present study we concentrate on the first source of scaling violation, namely, focus on 
the description of the perturbative evolution in the first two orders of the perturbation series in 
the QCD coupling at leading twist level, for the so-called skewed parton distributions (SPD) which 
arise in a number of processes such as deeply virtual Compton scattering [I], |5| , electroproduction 
of mesons @, || etc. The arguments in P, f|, ||, [7| convince us that these processes indeed factorize at 
this level into a perturbatively calculable hard scattering amplitude and a SPD. At twist two level, 
the dominant parton configuration in a hadronic wave function is that of two quarks or gluons. 
The main difference in the partonic description of these processes as compared to conventional 
deep inelastic scattering (DIS) is due to the fact that the longitudinal momentum fractions of 
particles propagating in the t-channel are different from one another and this difference is called 
skewedness r\. 

Unfortunately, a direct extraction of the SPDs from experimental data will hardly be possible 
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in general, since the processes are predicted as a convolution even at leading order (LO). Therefore, 
to deduce them one has to confront data with different models which are usually given at a much 
lower scale than factorization is expected to work and data are taken in experiments. Moreover, 
one of the central issues in the future should be to test the scaling violation phenomena in the 
cross section via evolution, similar to DIS. Thus, the main task is to describe evolution effects of 
SPDs as precise as possible. This includes the understanding of perturbative corrections at LO 
and beyond. 

The SPD is defined as a Fourier transform to the momentum fraction space of a light-ray 
operator constructed from yj-parton fields and sandwiched between hadronic states non-diagonal 
in momenta, schematically given by 



(x, V ,Q) = -L / dz.e ixz -{h(p')\0(z_,Q)\h(p)} with 0(z-,Q) = ^(0)^(z_ 

ZTT J 



Q 



Such a light-cone operator, 0(z-, Q), is a formal resummation of the usual local ones with definite 
twist. Its scale dependence is governed by a renormalization group equation (RGE), but with 
anomalous dimensions replaced by integral kernels depending on the positions on the light cone. 
Therefore, a SPD obeys an evolution equation with kernels given by integral transformations of the 
kernel in the light-cone position representation. However, the generalized skewed kinematics for 
the perturbative kernels can unambiguously be restored || from the conventional exclusive one, 
known as Efremov-Radyushkin-Brodsky-Lepage (ER-BL) region rj = 1. Throughout the paper we 
formally deal with ER-BL type equations || [10| 

d ., r 1 



dlnQ' 



-(j>(x, Q)=J q dyV (x, y\a a (Q)) 4>{y, Q), (2) 



where <fi = is the two-dimensional vector of the quark and gluon distribution amplitudes 
which mix under renormalization and V(x, y\a s ) for vector or axial-vector distribution amplitudes 
is a 2 x 2-matrix of evolution kernels given by a series in the coupling: 

00 \< 

V(s,y|a.) = 5:(-i) V®(x,y). (3) 

Until recently, only one-loop exclusive evolution kernels were available ||, [TD|, [IT], [T2|, |T3| and 
apart form the nonsinglet two-loop result [151 pi [T7| nothing has been known about higher order 
kernels. In our previous studies we have been able to calculate the two-loop anomalous dimensions 
matrices of the so-called conformal operators which are in one-to-one correspondence with the 
given distributions and finally obtain the singlet kernels in the form of exclusive convolutions. In 
the present study we conclude with giving, in great detail, the formalism we have used as well 
as presenting explicit results for the next-to-leading (NLO) kernels. As has been said above, a 
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simple continuation || allows one to easily deduce the generalized non-forward evolution kernels 
from the ER-BL ones. 

Different methods have been offered so far to solve the nonforward evolution equation. Two 



of them rely on direct numerical integration [rj| [19| of integro-differential equations. However, 
the first of them was designed for the treatment of the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi 
(DGLAP) region only, while the second one does not resum the InQ 2 terms. Next, one can 



reconstruct a SPD from its expansion w.r.t. an appropriate basis of polynomials ^T], with 
the expansion coefficients expressed in terms of conformal moments, which do not mix under 
renormalization in the one-loop approximation. This latter method can be employed to describe 
the evolution for general kinematics, however, it becomes inefficient in the small 77 region and in 
the t ~ 77 domain where SPDs rapidly change their shape. Note that the accuracy of this method is 
under control, however, the summation of a sufficient number of polynomials is time consuming. 
Another possibility is to use the conformal covariance in order to map the SPD to a forward 
distribution p3, which has in general no physical meaning, and then solve the corresponding 
DGLAP evolution equations with one of the standard methods. Yet another way is to solve 



the evolution equation in configuration space in terms of non-local conformal operators |25|, P6 



however, in this case, to evolve a given function explicitly the authors have used the same method 
of orthogonal polynomial reconstruction as used by us in previous studies. 

Beyond LO all methods except for the first two can be used in a straightforward manner only 
in an unrealistic conformal limit of QCD when the Gell-Man-Low /5-function is set equal to zero 
and by making use of the conformal subtraction (CS) scheme |27|, |28j which removes the special 
conformal symmetry breaking anomaly appearing in the minimal subtraction scheme, or they can 
be approximately applied for small skewedness, when the conformal non-covariant piece will die 
out. Thus, in order to be able to describe a realistic situation, we are left with a direct numerical 
integration of the evolution equations or with the orthogonal polynomial reconstruction method. 
Since all anomalous dimensions in NLO of conformal operators have been available for some time 



29], [H], |31|, we studied the evolution of several models of SPDs and have gained a first insight 
into the magnitude of evolution effects beyond LO for both the flavor non-singlet |32| and singlet 
sector [ 33j . However, due to complications as mentioned above, one should develop a more efficient 
numerical treatment which can be achieved by direct numerical integration routines. Therefore, 
the corresponding evolution kernels, whose Gegenbauer moments define the anomalous dimensions 
[p9| , p0| mentioned earlier, are needed in two-loop approximation. 

The analytical structure of these kernels is expected to be more complex than for the DGLAP 
ones. Therefore, a direct diagrammatical calculation in NLO is envisaged to be very cumbersome. 
Although, there is experience in the calculation of the ER-BL kernel in the flavor non-singlet 
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sector at two- loop order [|TJ, [16], |T7|], no appropriate technology has been developed so far for 
the remaining calculations. A technical analysis of the mathematical structure, which appears to 



be completely analogous in all channels, is given in [17], |3J]. Some particularly simple parts, i.e. 
renormalon chain contributions, of the non-singlet kernels were analyzed anew in Ref. ]35 . 



The goal of this paper is to present the method we employ for the reconstruction of the 
evolution kernels in NLO. It is based on the analysis of general properties of the kernels and 
known Gegenbauer moments. Assuming the CS scheme, which implies that the (local) conformal 
twist-two operators do not mix under renormalizationP], the anomalous dimensions are the same 
(up to a normalization factor) as in the forward case for DIS. This relation implies a one-to- 
one mapping between the evolution kernels of SPDs and the DGLAP-kernels. In the conformal 
limit of this special scheme, the reconstruction of the ER-BL kernels is reduced to an integral 



transformation from the known DGLAP-kernels [29]. The realistic situation is not as simple 



as that and beyond LO conformal covariance is broken in the conventional scheme, i.e. using 
dimensional regularization and (modified) minimal subtraction (MS), but it presents only a slight 
complication as compared to the restoration of the conformal covariant part of the kernels from 
the splitting functions. 

The outline of the paper is the following. We collect the properties of the ER-BL kernels in 
section ^ their support as well as relations to the DGLAP kernels, consequences of conformal 
constraints in QCD and conformal symmetry in M = 1 supersymmetric Yang-Mills (SYM) theory. 
We analyze the implications of these results for the structure of the flavor non-singlet ER-BL kernel 
in the two-loop approximation and extend it to all other channels. In section [| we reconstruct 
first the whole crossed-ladder diagram contributions by means of conformal and supersymmetric 
constraints from the known non-singlet results, then we restore the remaining diagonal pieces 
making use of the known splitting functions and present the results for the evolution kernels in 
the chiral even sector as convolutions of simple kernels having a one-loop structure. In section 
[| we perform the convolutions and present our results in an explicit form for the ER-BL type 
representation and extract from them the so-called skewed DGLAP kernels, which are needed for 
the numerical treatment of evolution in the DGLAP region. In section |5] we give a summary and 
conclusions. Different appendices are attached to give technical details. 



1 If we neglect non-covariant terms proportional to the Gell-Mann-Low /3-function which are induced by the 
trace anomaly, the conformal covariance can always be ensured by appropriate normalization conditions for the 
renormalized operators. 
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2 Structure of evolution kernels. 



In the following four subsections we give a survey of general properties of the twist-two singlet 



evolution kernels. Their support properties are pointed out in the next subsection 271 , which 
contains a generalization of the results obtained by Geyer, Robaschik, and collaborators in the 
past for the flavor non-singlet case || |36|, |37]]. Let us now give a few results which will be 



used later on: (i) The non-forward evolution kernels of SPDs are uniquely defined by an extension 
procedure from the ER-BL kernels, (ii) A simple limiting procedure reduces these generalized ER- 
BL kernels to the DGLAP kernels. In subsection |2.2| we demonstrate that if conformal symmetry 
holds true in perturbative QCD, the ER-BL kernels can be deduced from the DGLAP kernels 
by an integral transformation with a well-defined resolvent. We then discuss the breaking of 
conformal symmetry in the MS scheme and explain the third important issue, namely, (iii) how 
one part of the total evolution kernel is induced by special conformal anomalies. It is shown how 



to evaluate these anomalies in the ER-BL representation. In subsection |2.3| , we explain the use 
of the M = 1 SYM theory as a meaningful tool for the construction of the ER-BL kernels. Based 
on the results in the previous three subsections, we give a method for the construction of the 



two-loop ER-BL kernels in subsection [2.4| and demonstrate it in the flavor non-singlet case. The 
main result of this section is a representation of the NLO kernels in terms of convolutions of simple 
kernels possessing a one-loop structure, which turns out to be valid in all channels and allows one 
to restore the missing information from the known DGLAP kernels. 

2.1 Support properties of evolution kernels. 

In this subsection we recall the support properties of evolution kernels for SPDs. Let us start with 
the definition of these distributions at leading twist-two level in terms of light-ray operators that 
are defined for the vector (V) and axial-vector (A) case by: 



(: ° V \,.. .. s _ n . ,.. A 9* 



«i,« 2 ) = ^i(K 2 n) < \i)i{K X n), (4) 

7+75 I 



GQAj 



k u k 2 ) = G a + U 2 n){ \G a v+ {K 1 n), (5) 



where for brevity we have omitted a path ordered link factor that ensures gauge invariance. Here 
n and n* are light-like vectors with nn* = 1 and the plus and minus components of a four-vector 
are obtained by contraction with n and n*, respectively. Expansion of Eqs. ([|,[5]) into a series of 
local operators immediately guarantees the contribution from the leading twist-two only. They 
have even chirality and even (odd) parity for the (axial-)vector case. 

The SPDs are Fourier transforms w.r.t. the light-like distance between the fields of the light-ray 
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operators (f|,|5]), which are sandwiched between off-diagonal hadronic states: 

= jlj/ ^e^(P 2 5 2 | |^} («,-«) (6) 

given at a renormalization point [l. It can be shown that the support of the SPDs is \t\ < 1. Here 
the conjugate variable i plays the role of a momentum fraction and the skewedness parameter 
i] = A + / P + is defined as the +-component of the momentum transfer A = Pi — P 2 normalized to 
the +- component of P\ + P 2 , i.e. P + = n.(P\ + P%). More precisely, the momentum fractions of 
the incoming and outgoing partons are j^P+ and j^P + , respectively. Note that the SPDs also 
depend on the momentum transfer square A 2 , which, however, is irrelevant for the evolution and 
thus will be suppressed in what follows. A form factor decomposition would give us the functions 
H, E and H, E for the parity even and odd sectors, respectively, as introduced in ref. [0, which 
are governed by the same evolution equations as the distributions (§). However, we choose the 
normalization^ in such a way that in the forward limit A qf (t = z) they coincide for z > with 
the parton densities qf(z) and zg r (z) for A — Q, G, respectively, while z ) with z > are 

interpreted as antiquark distributions. 

The evolution equations for the SPDs (^J) arise from the renormalization group equations 
(RGEs) of the light-ray operators (|J) and (j5|). In the following we will discuss the flavor singlet 
case, where the quark and gluon operators mix with each other. For brevity we introduce the two 
dimensional vector 

and Nf is the number of active quark flavors. Note that due to Bose symmetry, the gluon operator 
also has definite symmetry with respect to the interchange of k% «-> k%, i.e. it is (anti)symmetric 
in the case of (axial-)vector couplings. The properties of (the non-local version of) the anomalous 
dimensions were intensively studied in the past on general grounds. For instance, the scaling and 
translation properties of the operators tell us that the general form of the RGE reads [37[] : 

k 2 ) = - J dy J dz~f(y, z; k 2 - fti) O («i[l -y] + K 2 y, k 2 [1 -z] + k x z) , (8) 

where the 2x2 matrix valued kernel 7 explicitly depends on the light-cone position k 2 — k% in 
the mixed channels: 

( z . ^ _ R s = ( Q °V(V> z ) ^ QG T(2/, z) \ (g) 

V ^ GQ j(y,z) GG i(y,z) )' 



2 For the gluon distribution we have included an extra pre-factor of 2 in comparison with the definition given in 
Ref. |. 
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The first important issue to understand is the anatomy of the evolution kernels for SPDs, which 
arise from the support property of the kernels 7. This problem can be solved by means of the 
a or Feynman-parameter representation of Green functions with a non-local operator insertion. 
It is sufficient to work in light-cone gauge and to formally generalize the a-representation for the 
gluon propagator JTJ. From these studies one can deduce the support of the kernels shown in Fig. 
1(a): 

7(2/, z; k 2 — «i) 7^ 0, for < y, z < 1; 0, otherwise. (10) 

Invariance under charge conjugation implies the following symmetry relation 

l{y, z; k 2 - Kx) = j(z, y; k 2 — «i). (11) 

It is also worth noting that the symmetry properties of the flavor singlet operators w.r.t. the 
interchange of their light cone arguments, i.e. K\ <-> k 2 , can be used to map the region y + z > 1 
into 1 > y + z by the substitution y — > 1 — z and z — > 1—y. Here the region 1 > y + z corresponds 
in the forward case to quark-quark mixing as it is the case in LO, while y + z > 1 appears due to 
a quark-antiquark interaction. 

From the definition of a SPD (|6]) and the RGE (H) we easily derive the evolution equation: 
d r 

H-i-q(t, r);n) = - dt'j(t, t', 77; a s (jj)) q{t', rj; p). (12) 
dfi J 

Here we use a definition of the singlet distributions, which is motivated by the one of normal 
parton densities, 

where the quark sector contains the sum of "parton" and "anti-parton" distributions. It is a 
simple exercise to see that the kernels for the SPDs are obtained by the integration 

7 (t, t 1 , rj) = t dy C dz ( QQl ^' Z) , QGliy ' Z)dt )s{t- t'(l - y - z) - (y - z) V ) , (14) 

where d t = d/dt and d^ 1 = f dt. Here the indefinite integration limits in the GQ channel 
induces an ambiguity which affects the unphysical moments only and has to be fixed by hand, 
e.g. by comparison of moments calculated in both representations. Note that this ambiguity 
implicitly appears also in the diagrammatic calculation of Feynman diagrams in the light-cone 
fraction representation and is responsible for different results given in the literature. It is worth 
mentioning that the representation fll4]) implies a simple scaling relation: 
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so that the entries are in fact two- variable functions of ratios. 

The invariance under charge conjugation implies now the symmetry for diagonal AAr y(t,t') = 
AA ^i(— t, — t') and off-diagonal A / B, AB j{t, t') = — AB ^(—t,—t') elements, while the connection 
between the parton-parton and parton-antiparton regions makes itself apparent in the substitution 
t' — > —t'. As explained, for the flavor singlet case it is enough to consider the region < 1 — y — z 
of the support < y, z < 1. Then the integral representation ([14]) implies the support shown in 
Fig. H (b), or formally 

7 (t, = &(t, t')f(t, o ± { * ^ Zl } > = W ~ + *) - W ~ 0^(-i - *), (16) 

where (— )+ stands for (off-) diagonal entries. The entries of the matrix f(t,t') are given by 

"/(MO = PWi/= 1+t ~l 1+t ' )m *= 1 ~'~i 1 ~' V |, (17) 



2 



with d t (d t x ) corresponding to (GQ). 




Figure 1: Support of the singlet anomalous dimensions 7 in light-cone position (a) and fraction (b) 
representation. In (b) we only show the support which arises from 1 — y — 2>0in the light-cone 
position representation. A second contribution that comes from the region 1 — y — z < can be 
formally obtained by t' — > — t', however, in the flavor singlet case it can be reduced to the first one 
by means of symmetry. Here the entries of / ±± = /(±t, ±i') are defined by Eqs. (0) and (|18D. 

A first glance at Fig. [II immediately gives the impression that the whole kernel can be obtained 
from the region \t\, \t'\ < 1. Indeed, it was proved in |jj that the continuation is unique (see 



S 



appendix 0). For practical purposes it is sufficient to replace the 9 structure: 

0(*-f)|t[,[*|<i-> (19) 
Thus, the evolution kernel for the corresponding SPD can be considered as a generalized ER-BL 
kernel and its restoration from a given ER-BL kernel is simple. 

There are further consequences arising from the evolution equation (|T2"P. If we replace in the 
definition @ for SPDs the off-diagonal hadronic states by diagonal ones, we immediately obtain 
the definition of the usual forward parton distributions (up to an additional 1/z for the gluon 
density) with t = z and r] = 0. Comparing the evolution equations, the DGLAP kernel appears 
as a limit of the generalized ER-BL kernel: 

i / QQy IQOv\ cxt / 1 \ 
PW = -7(^^ = 0) = -LIM 7 (f,O^-ton R (^| oc 7) (-,-). (20) 

An alternative derivation of this limit using moments is presented in appendix |B[ 

As we have just shown, it is sufficient for the following considerations to work in the exclusive 
kinematics. Thus, from now on we will deal only with the momentum fraction x — (1 + t)/2 and 
the common definition of ER-BL type kernels: 

V(x,y) = -~f(x-x,y -y)o<x, y <i, (21) 
here and everywhere else we use x = 1 — x. The general structure of the entries is 

AB v(x,y) = e(y-x) AB f(x,y)+9(y-x) AB g(x,y)±\ X ^ X _) for H * (22) 

As discussed above, the second ^-structure, i.e. 9(y — x), in the singlet case, can be removed by 
the (anti-) symmetry of the parton distributions w.r.t. y — > y. Again, the extension of the kernel 
in the whole region is done by a simple replacement of 9 functions, e.g. 

9(y - x) ^9(l-^) 9 (^\sign( y ). (23) 
2.2 Conformal properties of evolution kernels. 

To understand the classification of different contributions to the evolution kernel V with respect to 
the conformal transformation most clearly, we deal in this subsection with the so-called conformal 
operators J C?Jj and their anomalous dimension matrix 7. These operators build an infinite dimen- 
sional, irreducible representation of the collinear conformal algebra so(2, 1) in the space spanned 
by the bilinear field operators. This algebra arises from the full conformal algebra so (4, 2) by 
projection onto the light cone and consists of the generators of dilatation T>, special conformal 
transformation^ /C_, boost along the light-cone A4-+, and translation V+. Conformal operators 
3 The special conformal transformation is given by the product 1ZV C 7Z-, where the inversion 1Z acts as IZx^ = 
x^/x 2 and V c is a translation with a vector c. 
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are generated from the light-ray operators (|]) and (||) by differentiation w.r.t. k%, k,2'- 

^ = (id K1 + l d K2 y + ^ 2 c;lt\ A) - 3/2 (f^f;) ^ r (^ ^) lK1=K2=0 , (24) 

where v{Q) =3/2 and v(G) = 5/2. The index of the Gegenbauer polynomial CJ is determined 
by group theory: 2u(A) = 2d(A) + 2s(A) — 1, where d(A) and s(A) are the canonical dimension 
and the spin of the field of species A. The operators in Eq. ([24]) are the highest weight vectors 
/C_^C>J\ = and they carry the conformal spin j + 1 and the angular momentum / + 1. Acting 
with the generator of translation, i.e. applying on ^Jj, we generate the whole conformal 

tower of operators. In LO all members of a tower do not mix under renormalization, i.e. they 
have the same anomalous dimension 7j. This is a consequence of classical conformal symmetry 
and arises from the commutator constraints between the generators of dilatation T> and special 
conformal transformation /C_: [P,/C_]_ = i/C_, where as we have established above /C_ acts in 
the conformal tower step down operator. 

Scale symmetry is known to be broken by the trace anomaly in the energy momentum tensor 
[[38], [39], |4(| and is proportional to the Gell-Mann-Low /3-function. For the Green functions involv- 
ing only local field operators at different space-time points there is a one-to-one correspondence 
between breaking of special conformal and scale symmetries [[II] . However, as we will discuss in 
detail below, the renormalization of the Green functions with composite operator insertions there 
is an extra source of breaking of the special conformal covariance besides terms proportional to (3. 
Poincare invariance implies that the mixing matrix is triangular, thus, the RGE has the general 
form 

^=-£^°- M °*=X.U*,}- (25) 

To make contact with the (generalized) ER-BL kernel, we use the fact that the conformal moments 



of the SPDs (@) are given by the expectation values of conformal operators (24). For the flavor 
singlet case [see Eq. ([13|)] we have 



^ (A) - 3/2 £^^U-3 /2 (^ Y(t,m») = (P 2 ,s 2 \^o r jj \p 1 ,s 1 ),. (26) 

Moreover, since the Gegenbauer polynomials form a complete basis in the region [—1, 1], we can 
represent the kernel in the ER-BL region by the following sum 



'if 



j=v(A)-3/2 k=is(B)-3/2 JV i \ U \ A >) 



where w(x\u) = (xx) u is the weight function and Nj(u) = 2 Au+l '~fT^^f^r is the normal- 
ization coefficient. Up to an overall normalization, the conformal moments are identical with the 
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anomalous dimensions, i.e. AB Vj k = — AB jj k /2. Note that invariance under charge conjugation 
implies AB Vj k = for odd j — k. 

Fortunately, the additional special conformal symmetry breaking alluded to above is scheme 
dependent and as mentioned in the introduction it can be avoided in a special CS scheme. In 
such a CS scheme and in the formal conformal limit of QCD, when we set the Gell-Mann-Low 
function to zero, conformal operators do not mix with each other. However, once the /3-function 
is kept non-zero the RGE to all orders of perturbation theory is now 

= - § £ 2 Ajk ° H - (28) 

«M g k = Q 

Thus, in the CS scheme and up to a term proportional to the /5-function the ER-BL kernel is 
diagonal w.r.t. Gegenbauer polynomials: 

AB A^ v) = E - *) A % c-i%^ {B) {y - i/), (29) 

where the sum starts at j = for AB = QQ in the parity even case and j = 1 otherwise. The 
eigenvalues AB v l j are given by the diagonal entries, i.e. AB Vj = AB Vjj. From the representation (|29| ) 
it necessarily follows that the diagonal entries have the symmetry properties 



Taking into account the equalities among Gegenbauer polynomials with shifted index v, 

i^n^-^cr^ ( 31 ) 

we can easily derive necessary conditions for the kernels in the mixed channels, too: 

ox oy ox xx oy yy 

Since the matrix of conformal moments Vj are given by the forward anomalous dimensions 
known from DIS up to NLO, 

fM (33 ) 

we may reconstruct the kernels from these information with the help of Eq. (|29|). This represen- 
tation as an infinite sum can be converted into an integral by means of the generating function 

w (x\u) 



G(x,y;z\u) = E 1^TX C J( X ~ ^ C ^ y ~ & 
j=o iV i \r) 



T{v)T{v + l) 2 A "- l {xxy- l l 2 {l-z 2 ) 
TH^Tpu) [1 - 2 ((x - x)(y -y)- A^y) z + z 2 ] v+l 
, , , lQ^xxyyz 

X2-T1 



(34) 



u+l,u 
2v 



1 — 2 ((x — x)(y — y) — Ay/xxyy) z + z 2 
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which was obtained by Gegenbauer's summation theorem ||!J. For completeness we give here a 
simple reduction formula that is obtained from G(x,y; z\3/2) together with Eq. fl3lD : 

,1 ( QQy 19_QGy \ 

V D (x,y)= dz\ n 6d \ - \(z)G(x,y;z\3/2) 1 (35) 

\-\H x dx GG V -^$*idxj^ GG V ) y 1 y ' ' ! y 1 

where AB V (z) are related to the DGLAP kernels, however, defined as a Mellin transformation of 
conformal moments: 

1 rc+ioo 

V{z) = — djv jZ -i +1 . (36) 

ZTTt Jc— zoo 

Note that in the GQ and GG channel the first term of the expansion w.r.t. Gegenbauer polyno- 
mials does not exist and must be removed by hand. A different version of the reduction formula 
avoiding this problem has been applied at LO for the restoration of the ER-BL kernels from the 
DGLAP ones |]29|| . The evaluation of the integrals has been reduced by the Cauchy theorem to 
the evaluation of residues of simple poles, therefore, it was essential that the one-loop DGLAP 
kernels do not generate cuts in the complex plane. Beyond LO this property is lost and the derived 
reduction formula fl35|) is extremely hard to handle. 

To calculate the complete ER-BL kernel in any scheme, e.g. in the MS scheme with dimensional 
regularization, we have to take into account the breaking of conformal covariance induced by the 
conformal anomaly in the actionF[ There are different ways to account for this breaking: 

• An explicit calculation of the anomalous dimensions or kernels. However, depending on 
the representation, the calculation and/or extraction of the conformal non-covariant piece 
is extremely difficult or at least not straightforward. 

• Taking normalization conditions that ensure the conformal covariance of the renormalized 
operator w.r.t. the special conformal transformation for vanishing Gell-Man-Low function. 
For instance, doing so in LO, the anomalous dimensions are diagonal in NLO and the result 
in the MS scheme can be obtained by a finite renormalization extracted in LO. 

• One can analyze the conformal symmetry breaking with the help of conformal Ward identities 

(CWI) and derive constraints for the appearing anomalies which in turn will fix the piece 

we are interested in. This approach provides us also information on terms proportional to 

the Gell-Man-Low function. 

Beside the trace anomaly there also appear equation of motion terms and exact BRST operator insertions. 
In 4 — 2e dimensions the trace anomaly is proportional to (3 t (e,g) = —ge + (3(g) and we have to renormalize the 
product of trace anomaly and composite operators. Thus, this procedure together with the e dependent term 
causes anomalous contributions, e.g. anomalous dimensions in the case of dilatation. 
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All of the approaches sketched above are equivalent in the formal conformal limit. In order to 
reduce the effort in obtaining the desired result we prefer to work with the third one. Moreover, 
it enables us to get a deeper insight into the structure of conformal symmetry breaking countert- 
erms. It offers as well the technical tools to get the off-diagonal anomalous dimensions and the 
corresponding kernels in n th order by calculating Feynman graphs at (n — l) th order. For the 
reader's convenience we quickly outline this approach. The exact technical steps are published in 



great detail elsewhere |28|, |29], |30], |31], [43 . 

The major steps consist of: 
(i) Derivation of conformal Ward identities for the Green functions with conformal operator in- 
sertion 

Qjl{xx,x 2 , x N ) = (O j i(0)(f)(x 1 )(f)(x2) • • • <p{x N )) (37) 

by means of the path integral using dimensional regularization. The renormalization of the Ward 
identities provides us with a prescription for the calculation of the dilatation and the special 
conformal anomalies. Here we present these Ward identities in a simplified form required just to 
demonstrate the way conformal anomalies appear 
N ( d \ 

E \- X H^Z d *i Qjl{Xl,X 2 ,...,X N ) 



3=1 \ ® X Hj 



.1 



J2 {(I + 3)1 + 7} jk Qki(xi,x 2 , ...,x N ) + ..., (3£ 



k=0 



N ( d d \ 

H U l y 2X M Xv i]fa~. + X ^dx~. ~ 2 ( d +* + S ^ X ^J ^( X 1' X 2, ...,X N ) 

= £ {a(j, 1)1 + 7 c (0} jfc 3ki(xi, x 2 ,...,x N ) + ..., (39) 

where the anomalous dimensions matrix 7 fe and the so-called special conformal anomaly matrix 
7j fc are induced by quantum fluctuations. We introduced the conventions a(j, I) = 2(j—l)(j+l+3) 
for the coefficient which appears at tree level in the transformation of a conformal operator under 
special conformal transformations. The dots in Eqs. (|38| , |3"9"D stand for terms involving Green 
functions with the renormalized product of a conformal operator and a conformal variation of the 
action ([O j i5S](f)(x 1 )(f)(x2) ■ ■ .(J)(xn))- 

(ii) As a next step we derive matrix constraints for the conformal anomalies and find their so- 
lution. It turns out that there are two commutator relations stemming from the algebra of the 
collinear conformal group which provide non-trivial relations between the conformal anomalies. 
For vanishing Gell-Mann-Low function (complete constraints) we have 

[P,/C_]_ = */C_ =}► [a(/)+ 7 c (0,7L=0, (40) 

[V + ,IC_]_ = 2z(V + M- + ) =► 7 c (^ + l)-7 c (0 = -27- (41) 
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Since the matrix a is diagonal, a,jk = a(j, l)5jk, the constraint (f40|) tells us that the off-diagonal ma- 
trix elements of the anomalous dimension 7 are indeed induced by the special conformal anomaly 
matrix 7 C (/), while the spin dependence of 7 C (/), induced by the breaking of Poincare invariance 
due to the special conformal transformation, is governed by Eq. fl4"ip. Note that the diagonal 
form of the special conformal anomaly 7 C (/) is necessary and sufficient for the vanishing of the 
off-diagonal part of the anomalous dimension matrix 7. This ensures the existence of the CS 
scheme, in which the conformal covariance of the operators holds true in the conformal limit, 
(3 = 0, at any order of perturbation theory [pTjl - 

The extension of these constraints ( f40| ) to nonvanishing /3-function is straightforward, however, 
they require additional algebra [E5[ [ft], [JTJ] , which results in the following change of Eq. 



a(Z)+ 7 c (Z) +2^>(Z),7 
g 



0. 



(42) 



whereQ b jk (l) = 6 jk [2(1 + k + 3)6 



[1 + f- 



1) J k ](2k + 3) j 1. Decomposing the anomalous di- 



mension matrix in its diagonal (D) and off-diagonal (ND) part, the solution of this constraint can 
be constructed by successive approximations: 

7 c (0 + 2fb(/),A 



7 ND (<?) 



7 D (S) = -<?7 D Q?) + ---, with QA 



jk 



a(j, k) 



(43) 



(iii) The last step consists of the explicit evaluation of the anomalies. The LO anomalous dimen- 
sions of conformal operators have been known for a long time [T2|, 1221, |2B 



QQJ 3 ) 



'j 



GQJ°) 
13 



-cv 



(j + l)(j + 2) 



-24N f T F 



j'0' + l)0' + 2)0' + 3) 



-4V0' + 2)+4V(l) 

j 2 + 3j+4, 
j(j + 3), 



3(j + 1)(j + 2) 



J 2 + 3j+4, 

i(i + 3), 



-C A -4V(j + 2)+4V(l)- 



8C A 



j(j + l)(j + 2)(j + 3) 



x 



J 2 + 3j + 3, 
. j(j + 3), 



where the upper (lower) row corresponds to even (odd) parity and (3q = ^TfNp — 
special conformal anomalies we found at one-loop level |J2^, |3y, |31[ : 

7 c(0) = _ b7 (0) + w 
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We use the following definition of the step-function: 6jk = {1 5 f° r j — k > 0; 0, for j — k < 0}. 



(44) 
(45) 
(46) 
(47) 



C A . For the 



(4J 
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The first term on the r.h.s. appears as a result of the one-loop renormalization of a conformal 
operator in a subgraph and we observe that it is induced by the breaking of scale invariance. For 
the Green functions constructed from local field operators this exhausts the sources of symmetry 
breaking, however, for the case at hand it is no longer true and we have an addendum w. The w 
matrix contains new information from the renormalization of the operator product of conformal 
operators and the conformal variation of the action in 4 — 2e dimensions. These renormalization 
constants are calculated by means of modified Feynman rules. It turns out that the w matrix is 
universal for vector and axial-vector operators and its matrix elements read 



9^ 2<k (3 + 2k) 



QG„ 



= -2c F [l + {-iy- k 

x l2A jk + (A jh - VO* + 2) + V(l)) 
= 



GQ 



-2C f 



l + (-iy- fc R-_ 2 , fc (3 + 2A;) 



(j-k)(j + k + 3) \ 
(k + l)(k + 2) J ' 

l(j-k)(j + k + 3) 



GQ 



Wjk 



= -2C A [l + (-iy- k ]9^ k (3 + 2k) 
x \2A jk + (A jk -ijj(j + 2)+7p(l)) 



6 (* + !)(*; + 2) 



T(j + 4)r(fc) 

T(j)T(k + A) 

The elements of the matrix A are rather complicated 



(49) 
(50) 
(51) 
(52) 



+ 2{j-k)(j + k + 3) 



T(k + A) 



A jk = v ( 3 + 1 l +A ) - v» ( J -^) + m -k)- v»o" + 2) - v(i). 



(53) 



It is instructive to note that the special conformal anomalies (fHf) obey certain constraints which 
originate from the anomalous superconformal WI and a commutator of the generators of the 
special conformal transformation and restricted supersymmetry, and imply that one can restore 



all entries in Eqs. (|4"9"D- (|52"D from the knowledge of e.g. ^7 and superconformal anomalies [M\ 

The form of the special conformal anomaly ( 
diagonal part of the anomalous dimension matrix 



implies the following structure for the off- 



7 



ND(0) 



7 (0) ,(A)l-7 (0) )rf + <7 



(54) 



where we defined the matrices djk = bj k /a(j,k) and g^ k = Wjk/a(j, k). This is the final result 
of the conformal approach which give together with the known anomalous dimensions from DIS, 
see for instance |45| , fj6| , the complete anomalous dimension matrix to NLO in the MS scheme 
relevant to non-forward processes. 

As was explained in the introduction it is paramount to know the ER-BL kernels. To obtain 
them, one can transform, in a straightforward way, the conformal anomalies into the momentum 
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fraction representation. For instance, for the QQ channel this kernel reads 



\ QQ w(x,y)] , = QQ w(x,y) - 8(x - y) / dz QQ w(z,y) + —5(x - y) / dz(z - y) QQ w(z,y), 

Jo dx Jo 



QQ w(x,y) 



-Cp9(y — x 



x 



y(x- yf 



x 
V 



x 

y 



(55) 



In the next step we construct the ®®g(x,y) kernel with its conformal moments satisfying the 
equality 2[(J + l)(j + 2) — (k + l)(k + 2)] QQ g jk = QQ Wjk- Using the eigenvalue equation for 
Gegenbauer polynomials, 



dx 2 



xxCf 2 (x - x)] = -{j + l)(j + 2)C| /2 (x - x), 



(56) 



and the representation (|27|), we can easily write down a second order differential equation for the 
®Qg(x, y) kernel 



a 



(J 2 



xx— QQ g(x,y)-— yy QQ g(x,y) 



dx 



dy 2 



\)(x, y), for x 7^ y. 



(57) 



The solution of the homogeneous equation is purely diagonal w.r.t. Gegenbauer polynomials, thus, 
we find from Eq. ( |57|) the ®®g kernel, which contains beside the desired off-diagonal also (arbitrary) 
diagonal conformal moments. Similarly one constructs the other channels. Our results for the 
whole singlet sector reads 



g(x,y) = 6{y - x) 



( 



V 



-cv 



In 1- 



v I 



y-x 



-c. 





lnf 1- 



y—x 



\ 



J 



± 



X — > X 



(58) 



with the (— )+ sign corresponding to (non-) diagonal entries and the "-(-"-prescription defined as 
[V(x, y)] + = V(x, y) — S(x — y) Jq dzV(z, y) + const • S(x — y), where the constant term is fixed in 
appendix |B[ 

Now we come to the dotted kernels v(x,y) which posses the conformal moments y(°\d . 
The matrix djk can be generated by a derivative w.r.t. the index of Gegenbauer polynomials 



d q(t)i,= 3/2 = -2 E d jk cl /2 {t), ^c]_ x (t\ v= , l2 = -2 £ d Jk cl\(t), 

k=o av k=l 



du 



(59) 



and analogous relations, but with a different sign and range of summation on the r.h.s., follow 
immediately from the orthogonality relation, for the function w {x\i>) Cj(x — x)/Nj (u). Thus, it 
is obvious that the dotted kernels can be generated from 



w {x\u(A) + e) c v( A )+e 



AB v*(x,y\e) - N _ {v{A) + e) ^ +3/2 - uiA) 



ABi(0) r v(B)+e 

03 °j+3/2-i/(B) 



{x-x) AB vfC 



(y-y) 



(60) 
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by differentiation w.r.t. the parameter e at e = 0. Of course, for diagonal entries the generating 
kernel AA v z (x, y\0) is symmetric w.r.t. the weight function AA w(x\v) [see relation (|30|)1. Thus Eq. 
d60|) provides essentially a logarithmic modification of the LO kernel 

A %\x,y)=6(y-x) A Y(x,y)±l X y 2l} for {a^B' (61) 



However, using the symmetry relation (|32|) to get the mixed channels one has to fix the integration 



constant appropriately. Thus, the generic form of the dotted kernel is 

AB v\x,y) = 6(y-x) A Y(x,y)\n*+A AB r(x,y)±{ X ^*\ for H^. (62) 

y [y ^y J [ A r B 

To obtain the addendum A AB f\ it is necessary to have a closer look at the structure of LO kernels, 
which will be given in the next subsection. 

The factorized structure of NLO off-diagonal anomalous dimensions is transferred to the mo- 
mentum fraction kernels and can be constructed out of the conformal anomalies we have just 
found. In convolution form they read 



V ND W(x, y) = -{X-V){v®( V® + + 



2 



9 



(x,y), (63) 



f 1 

where ri®r 2 (a;, y) = / dzrAx, z)T 2 (z,y) defines the exclusive convolution and [X — V) projects 
Jo 

out the diagonal part. 

2.3 jV = 1 supersymmetric constraints. 



Since the two-loop off-diagonal kernels fl6"B| ) are available we are left with the construction of their 
diagonal counterparts. Unfortunately, the direct use of the integral transformation (|35| ) with the 
generating function given in Eq. ([34]) is difficult due to the complicated analytical structure of 
the integrand in the complex plane: cuts appear on top of poles. Therefore, we have to look for 
complimentary sources of restrictions on the form of the NLO kernels. Sufficient information comes 
from the constraints deduced from the graded commutator algebra of J\f = 1 SYM theory |14|, [50 



There is a one-to-one correspondence between QCD and SYM Lagrangians provided one identifies 
quarks with gluinos and put the former in the adjoint representation of the colour group. We can 
map thus the QCD result to SYM theory by equating the colour factors: Cp = 2Tp = Ca = N c . 
The main use of the previously mentioned constraints will be in the reconstruction of the crossed 
ladder diagram contributions which have the most complicated analytical structure. For this type 
of Feynman graphs the restoration of colour factors is unique. 

Since the conformal symmetry breaking part has been previously fixed, we can legitimately 
assume in what follows that conformal covariance holds for the anomalous dimensions. In addition 
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to this, we have to use the renormalization procedure preserving supersymmetry. In reality none of 
these statements hold true in the MS scheme beyond LO. From the commutator of the dilatation 
and restricted supersymmetry generators [Q, X>]_ = ~Q applied to the Green function (|37|) one 
finds six constraints for eight anomalous dimensions in the chiral even case 

Q % + - G % = i Q % + G %, ^ = V,A (64) 



j Q 

QQ l] + i + ^ GQ l] + i = QQ lt~ 3 f G ^' and V&A (65) 



Hj ~ = 0, (66) 

J 6 J o 



which we call type-I (or Dokshitzer relation [J49|), II and III constraints, respectively. It is worth 
to mention that the first two constraints are valid to all orders of perturbation theory in the super- 
symmetry preserving regularization and renormalization scheme while the last one gets modified 
beyond one- loop order by conformal non-covariant, off-diagonal elements of the anomalous dimen- 
sion matrix |50f . 

As an example of the power of these constraints we demonstrate their use at LO. We show 
that all entries in the singlet sector can be reconstructed from the QQ channels by solving the 
constraints (|6^ ) -(l66| ) . At first we consider the anomalous dimensions, and eliminate with the help 
of the type-Ill relation the QG entries in the type-II relation. By elimination of the GQ entry 
either in the vector or axial- vector channel we obtain a recurrence relation for the remaining mixed 
channel. For instance, the solution in the parity odd sector reads 

G V = (67) 

^ (j + l)(j + 2) {0i> 



i=2,4,. 



where in addition the (analytic continuation to the non-physical region of the) lowest moment has 
to be known, i.e. G( ^)q- The remaining five entries are obviously linear combinations of the known 
ones: The type-II relation gives the QG entry for even parity, the typ-III constraints provide the 
remaining two mixed channel entries, and the GG channel follows from the Dokshitzer relation. 
In this way we obtain the whole anomalous dimension matrix for SYM theory in LO. Let us 
add for completeness that the restoration of colour factors is unique up to the GG channels, 
where the self-energy contribution provides a constant term proportional to Nf. This constant 
can be fixed either by explicit calculation of self-energy insertions in QCD or by the fact that the 
energy-momentum tensor is conserved and its anomalous dimension vanishes. 

The remaining goal is to derive and solve constraints for the ER-BL kernels by means of the 
representation ( p9|) . Unfortunately, we can only write down the type-II constraints in terms of 
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integral transformations such as given in Eq. fl35|). Therefore, we expect a cumbersome solution 
of these constraints, which is maybe useless for practical purposes. Fortunately, the knowledge of 
the eigenvalues give us a hint to overcome this problem as described below. The representation 
( |29"D together with the type-Ill relation (|66| ) provides us with the relation between the kernels in 
the mixed channels 

G % i (x,y) = ^- C ^) i (y,x). (68) 

yy 

Next we can eliminate the GQ entry in the Dokshitzer relation and take into account the relations 
fl3T|) between the Gegenbauer polynomials, to find a differential equation for the GG kernels: 

— QG v l (x,y) + — GG v i (x,y) = -3 QG v i (x,y) for i = V, A. (69) 
ay ox 

Thus, it remains to find equations which determine the mixed channel contribution through the 
knowledge of the quark one. A closer look at the anomalous dimensions (^)-(^) shows that the 
QG entry for parity even and the difference between parity even and odd in the GQ channel, i.e. 
G Qry 5 = GQ^ V _gq^ v ^ are p r0 p 0r tional l/y[j+l)(J+2)j , which coincide with the eigenvalues of the 
scalar 0^ D=6 ) theory kernel G[G v a (x, y) = 8(y — x)x/y + 9{x — y)x/y, which appears as a part of the 
whole QQ kernel. Thus, we can write down for the kernels ® G v A as well as for ® G v 5 = OG v v — Q°v A 
the following differential equations 

Q G V A (x,y) = ^% a (x,y), (70) 

dy 

d d 

— GQ v 5 (x,y) = -4 QQ v a (x,y) + const • xx and — GQ v s (x,y) = 2 GG v a (x } y), (71) 



where the kernel GG v a (x,y) has the eigenvalues 2/^(j + l)(j + 2) J and can be constructed from 
the QQ ones by means of the differential equation ( |B"9"D with the r.h.s. set equal to zero. In the 
first equation of the set (|7ID we included a term proportional to xx, which reflect the fact that 
the conformal expansion of GG v s (x,y) starts with j = 1 [compare with Eq. Consequently, 
the lowest conformal moment of the r.h.s. has to vanish. Thus, in our case this constant is 
const = A Q %/N (3/2) = 12. 

Let us now discuss the construction of the LO kernels from the knowledge of the QQ entry 



QQ v(x,y)] with QQ V = QQ V " + QQ V \ Q%\ x , y) = 6(y - x)-—— + \ X ^ X . (72) 
J+ yy-x [y-*y\ 



In the parity odd sector Eq. (|7Dp give us the QG entry and Eq. (^) the GQ one. The integration 
of the constraint (69) provides us then with the missing GG kernel. The integration constant c(y) 



as a function of y is almost determined by the necessary condition (|30D for the kernels to have 
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a diagonal form. The remaining degree of freedom, i.e. {yy) 2 c'(y) = (xx) 2 c'(x) = const, can be 
easily fixed by the requirement that GG Vj 1 = for j > 0. Now we consider the difference between 
parity even and odd cases. The set of differential equations (ffT|) gives us the corresponding GQ 
entry, while the remaining integration constant is fixed from the requirement that the conformal 
moments G< ^ vanish for j > 0. Then the difference in the QG channel follows from the symmetry 
relation (j68|) . Finally, the integration of Eq. ( |69"D gives us the contribution in the GG channel. 

The results obtained, with a minimal calculation of Feynman diagrams, has the advantage that 
the kernels are diagonal for all conformal moments. A direct calculation suffers in the parity even 
case from subtleties, which generate in the unphysical sector off-diagonal conformal moments. 
Now we present the improved kernels in such a way, that the underlying symmetries are explicitly 
manifest. For / = {A, V} the kernels read 

, nl , / C F \ QQ u(x,y)} 2T F N f QG v I (x,y) \ 

\C F GQ v\x,y) C A GG v I (x,y)-f5(x-y) /' 

where the entries for parity odd are given by only two types of kernels and the difference to the 
even case arise from a third one: 

QQ V = QQyO- _|_ QQ y 6 ) QG y A _ _QG y a^ GQ y A _ GQ y a ^ GG y A _ ^ GG y a + GG^bl ^ /y^N 
Q%V = QG v A_ 2 QG v c^ GB v V ^GB v A + 2 GB v c for R = {g> G} + (?5) 

The functions AB v l are defined by Eq. fl6"T|) with 

j AB f a \ x*W*f 1 1 fA,B = {Q,G} 

AA,c ^ W ~ 1/2 [2gy[j-ln(ay)]+y-x1 fQ 

The index z/(v4) coincides again with the index of the Gegenbauer polynomials. It is important to 
note that we have introduced different "-(-"-definitions in the QQ and GG channels in order to have 
a one-to-one correspondence between exclusive and inclusive kernels which require a regularization. 
For the exclusive QQ kernels we use the conventional prescription 

[V(x,y)] + = V(x,y)-5(x-y) [' dzV(z,y), (79) 

Jo 

while for the GG kernel an extra term const (V)5(x — y) is subtracted, with the constant expressed 
in terms of an integral of the kernel itself [see Eq. (|206|) 1. 
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'fay) = - QG v c (x,y) 



(80) 



We also introduced in Eq. ( |T?D the ®% c kernel, which does not show up in the LO kernel, but 
is of importance beyond one-loop approximation. It satisfies the differential equation 

d_ 

dy 

and is diagonal w.r.t. Gegenbauer polynomials. The given explicit representation of this kernel 
will be essential for the analysis of the ER-BL kernel in NLO. 

It is worth to mention that the eigenvalues of the same iZ-kernel in different channels are 
related to each other (here Vjj = i>„- 

1, 
6 



QQ„,a 
3 



■ QG a 
3 



,j by the relations 

6 nn „ 1 



GQ a 

J(J + 3) j 



_GG a 

2 j 



GG b 
3 



2V>0" + 2 ) + 2^(1) + 2 

6 nr. „ 1 



LQG c _ 
6 ' j(j + 3) 



GQ v c 
3 



(i + i)(i + 2) : 

2 



- GG V C = 

3 1 i(i + i)0' + 2)(i + 3) : 



(81) 
(82) 
(83) 



where for the AA Vj kernel the proper "-(-"-prescription has been taken into account. 

Even though we do not consider the chiral odd sector in our presentation we just want to 
mention that the corresponding constraint can be easily implemented for the kernels 



^% T (x,y) + ^ GG v T (x,y) 



0. 



It allows us to find the GG kernel in a unique way from the QQ entry as discussed above: 



AA V i0)T (x,y) 



C F 
C A 



iA u\x,y) 



2 S ( x -y) 



c F 

4C A + (3o 



for A 



Q 
G. 



34) 



35) 



2.4 Method of reconstruction of two-loop kernels. 



In the previous sections we have established the complete structure of the ER-BL kernels in LO 
and also partially in NLO. For the MS scheme the piece containing the off-diagonal part to NLO 
is known as a convolution (|63D , however, we have no explicit representation of the projection 
operator (X — T>) at hand which makes the following analysis more complicated. Let us drop it 
and represent the whole contribution to NLO as 



V (1 \x,y) = -V 



fay) 



9 



fay) + ' D (x,y), 



36) 



where T>(x,y) is a pure diagonal part. We adopt the following strategy for our considerations. 
To restore this part one goes first to the forward limit (|2"DD in which all non-diagonal terms die 
out and compares the result with the known two-loop DGLAP kernels P^\z) in order to get the 
DGLAP representation of 72(x,y): 

T>(z) = UMD(x,y) (87) 



P {1 \z) -UM{-V- 



'y-(o) + A ^ 



fay) 



9 



y(°) 



fay) 
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Afterwards one may try to use the reduction formula fl35| ) to restore the ER-BL representation. 
Unfortunately, as we mentioned several times, this last step is too cumbersome to be performed in 
an analytical manner. Thus, we are forced to apply supersymmetry to get a better understanding 
of the structure and in order to reconstruct the missing diagonal pieces. If one looks to the DGLAP 
kernels, one immediately sees that the most complicated terms are Spence functions multiplied 
by the kernels appearing in LO 

x 



±2 AB p{-z)S 2 {-z)+ AB p{z 



ln 2 z 



2C(2) , where S 2 (z) 



V(i+«) dx , 1 
— In — 

z/(l+z) X X 



Such contributions arise from the crossed ladder Feynman diagrams. Since they have no ultraviolet 
divergent subgraphs and, thus, require no subtraction, it appears just as in LO only as a 1/e 
divergence. Therefore, conformal covariance and supersymmetric constraints hold true for these 
contributions. Of course, it is questionable if these six constraints relate all contributions of this 
set of Feynman graphs in a given gauge, e.g. light-cone gauge. Fortunately, the answer is irrelevant 
for the reconstruction of the most complicated part containing the Spence functions. Since the 



ER-BL kernel is known in the non-singlet case [15, 16, 17|, we have the one-to-one correspondence 
of this piece in the DGLAP and ER-BL representation for the QQ channel at hand |I| which helps 
in the reconstruction of the other channels. 

It is convenient to introduce the following notation for the NLO result, 



V^(x,y) = C F [v{x,y)\ 
V^(x,y) = C F 



+ ■ 



v(x,y) = 9(y - x)f(x,y) + 



x 
V 



x 

y. 



f 



2y a + QQj 



C F V F (x,y) - —Vp(x,y) 



C F -^)V G {x,y) 



(89) 



where the functions ®Qf with i = {a, b} are given in Eq. (ff6|). The piece which mainly originates 
from the crossed ladder diagram is proportional to (C F — C^/2), 



V G {x,y) = 2v a {x,y) + -v(x,y) + f G(x,y) + 



(90) 



x — > x I \ 

y^y)J ' 

and is diagonal w.r.t. the Gegenbauer polynomials. This is obvious for the first two terms appearing 
in Eq. (p0|). The function^ G(x, y) contains both ^-structures 

G{x, y) = 6{y - x)H{x, y) + 8{y - x)H(x, y), (91) 

where the functions H(x,y) and H(x,y) are build from the LO function f(x,y) and f(x,y) = 
f(x, y) in combination with the Spence function Li 2 and double logs: 



H(x,y) 
H(x,y) 



f(Li 2 (x) + \ny\nx) - fhi 2 {y) 



x 



(f - f) Li 2 1 - - + -hSy + / (U 2 (y) - U 2 (x) - \ny\nx) 



3 We have slightly changed the original definition given in |T^, |l6| by G(x, y) + 26(y — x)f\nylnx 



(92) 
, (93) 
G(x,y). 
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where Li2 (2/) = — Jq dx [hi(l — x)]/x. It can be easily checked that the G-contribution is symmetric 
w.r.t. the weight xx. Note that the two different ^-structures do not obey separately this symmetry 
and contain therefore off-diagonal conformal moments. At this stage it is not so important to 
redefine H and H so that they are separately symmetric and diagonal. This will be done for the 
presentation of the explicit NLO result in section |j. Performing the limit (|20| ) [of course, only 
the QQ channel is relevant] we obtain the following correspondence with the non-singlet DGLAP 
kernel 

G(z) = UMG(x, y) = 6(z)6(l - z)H(z) + 9(-z)9(l + z)H(z), (94) 

where 

H{z) = UMH(x,y)=p(z)(\n 2 z-2((2))+T(z), (95) 
H(z) = UMH(x,y) = 2p(z)S 2 (-z) +T(-z), T(z) = 2(1 + z)\nz + 4(1 - z). (96) 

Up to the simpler term T(z), we recover by mapping of — 1 < z < into the < z < 1 region the 
desired expression (|88| ). 

As demonstrated in subsection |2.3| for the LO kernels, we can now apply six supersymmetric 
constraints to find the G kernels in all other channels from the kernel (pl|) in the QQ channel. 
This task will be achieved in the next sections. 

Now we study the structure of the remaining diagonal piece D(x,y) in the non-singlet sector. 
Knowing the correspondence of the G kernels in different representations, the remaining diagonal 



piece of Vg(x } y) can be easily restored from the DGLAP kernel |5T 



P G (z) = 2p a (z) + ~p(z) + G(z), (97) 

where p a (z) = LlMv a (x, y) = 1 — z and p(z) = LlMv(x, y) = (1 + z 2 )/(l — z) by substitution 

p a (z) -> v a (x, y), p(z) -> v(x, y). (98) 

Next, the Feynman diagrams containing vertex and self-energy corrections provide Vp pro- 
portional to (3q. Its off-diagonal part is induced by the renormalization of the coupling and is 
contained in the dotted kernel: 

V p (x,y)=v(x,y)+D l3 {x,y). (99) 
This dotted kernel is given in Eq. (|6^ ), where Af is zero in the QQ sector: 

v(x, y) = 6(y - x)f(x, y)\n- + T ^ X ) . (100) 

y {y^yj 
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This can be proved easily by forming conformal moments and using the triangular form of this 
matrix. The remaining diagonal piece, Dp, is deduced from the known NLO DGLAP kernel [51 



5 

Pp(z) = -p(z) +p a (z) +p(z), p(z)=UMv(x,y)=p(z)\nz + l-z, (101) 

by going to the forward kinematics and restoring the missing contributions from it by the substi- 
tutions (H): 

Dp(x,y) = ^v(x,y) + v a (x,y). (102) 

Indeed, the final result coincides with ]T5|, [IB], [T7] . 

Making use of the known non-diagonal part in the flavor non-singlet sector [which is the same 
as appearing in the QQ channel of Eq. fl36f)1, V F can be represented up to a pure diagonal term, 
denoted as D F (x,y), by the convolution^ 

(e e e \ 

v®v + g®v - v®g J (x, y) + D F (x, y), (103) 

where the g kernel is given by the QQ entry of the matrix (|58|). To find an appropriate represen- 
tation of this missing diagonal element we first take the forward limit. Since the forward limit of 
the convolution is 

LIM hu]+ ®[v}+\ = {LIM [v]+} ® {LIM [v]+} , (104) 
where we have introduced the inclusive convolution 



P 1 (z)®P 2 (z) = / dx dy5(z-xy)P 1 (x)P 2 (y), 
Jo Jo 

with the commutator g <g> — ® g droping out in the forward limit, we obtain 

HMV F (x,y) = -pk>p + UMD F (x 1 y). (105) 
The comparison of LIMVp(x, y) with the corresponding part of the DGLAP kernel J5T]] 

P F (z) = ^-2((2)-^lnz + ln 2 z-2\nz\n(l-z)^p(z) 

1 — 3z , 1 + z 9 . 
+ l-z + — - — \nz — In 2 z, (106) 

yields the result in which all double log terms are contained in the convolution 



p®p 



27 3 
-— + 2C(2) + - In z - In 2 z + 2 In z ln(l - z) 

+ 2(1 - z) In 1 \n 2 z 



p{z) (107) 



1-z 



7 We remind the reader that all flavor non-singlet kernels are supplemented by a "+" -prescription. Therefore, 
we use for simplicity in this section the convention that C — A® B is indeed defined by the convolution of 
[C}+ = [A]+®[B]+. 
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and, therefore, only single logs survive in D F (z) = LlMD F (x,y): 

D F (z) = P F +pk>p 

= -\p\-z) \nz-p a (z) jlnz - 21n(l - z) - ^} - ^p(z). 



{10t 



Here we have introduced for convenience the kernel p a (z) = 1 — z. The next important point is 
that the remaining log terms can be represented as convolutions of p a and p, i.e. 



p a ®p a 



[3p a (z) +p a {-z) Inz] 



+ ' 



p ®p 
Thus, we finally have 



-p{z)+p a {z) 



hxz - 21n(l - z) 



D F {z) = ^p a ®{2p + p a } (z) + ^p(z) + ^p a {z). 



(109) 
(110) 

(111) 



Since D F (x,y) is by definition diagonal, the extension of D F (z) towards the ER-BL kinematics is 
trivial: 



D F {z) -»■ D F {x,y) = ^v a ®(2v + v a ) {x, y) + ^v{x, y) + ^v a (x, y). 



'112) 



Evaluating the convolutions one can establish the equivalence of this equation with the explicit 



calculated expression in Ref. [|T|, |T(| [Tfl| : 

f /4 

V F (x,y) = 6(y-x)< 



2C(2) / + 3 



x 



3 x 



x 



X 



In- -(/-/) In -In 1- 



x 



y 



X 



X 



+ / + — In - ) -\nx (1 + lnx - 21nx) + 



2y 



2y 



x — ► X 

y^y 



;ii3) 



Recapitulating our results presented in this section, the ER-BL kernel to NLO has a rather 
simple structure governed by conformal anomalies, the crossed ladder contribution G(x,y), as 
well as a remaining diagonal piece D(x,y): 



VW(x,y) = -V I ®lvW + ^lUx,y) 



0o 



9 



V® 1 



(x,y) + G I (x,y) + D I (x,y). (114) 



The first two terms on the r.h.s. are known, the kernel G(x, y) can be restored from the QQ entry 
by means of supersymmetric constraints, and the missing diagonal piece can be extracted in the 
forward case from the DGLAP kernels P^ 1 \z): 



DHz) = P {1)I (z)-UM 



V® 1 + —1 



9 



+ G 1 . 



115) 



In the non-singlet case we saw that this piece does not contain double logs or Spence functions. 
Moreover, it can be represented in a rather straightforward manner by simple kernels, known 
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from the LO result, and their convolutions. This immediately allows one to obtain the ER-BL 
representation. Since the established structure is deeply related to the topology of Feynman graphs 
and the renormalization of the subgraphs, we expect that the missing entry D\x,y) in all other 
cases can be build by means of known kernels, too. 



3 Reconstruction of evolution kernels in NLO. 

In this section we present, in detail, the reconstruction of all eight evolution kernels appearing in 
the vector and axial-vector case. In the following two subsections we complete the construction 
of V and of G 1 kernels, respectively. Afterwards we extract the remaining diagonal part and 
represent the ER-BL type kernels as convolution of kernels, appearing in LO or dotted kernels, 
and the G kernels. 



3.1 Construction of V. 

In this subsection we present a systematic construction of the dotted kernels **9y*(:r, y), which 
posses the off-diagonal conformal moments [ AB v] - AB v\ )djk- Their generic form is given in Eq. 



(|62"D and the remaining goal is to determine the addenda A AB f l (x,y). As we already saw in 
the previous sections, in the QQ channel there is no such extra term for the a and b kernels. 
From the generating kernel (^) for the QQ channel we can derive all the other ones. Since the 
eigenvalues of the a-kernel in the QG channel are ^® G Vj = —^°Vj, we can obtain the dotted kernel 
by differentiation w.r.t. y: 

QG v a (x,y) = —^- QQ u a (x,y) + ---, (116) 
dy 

where the ellipsis denotes possible diagonal terms. One can easily check that this equation is 
satisfied by Eq. fl6"2|) with A° :G f a (x,y) = 0. To find the dotted a-kernel in the GQ entry, we use a 
symmetry relation for the generating kernels, which arise from - QG v^/Q = Q GQ upj(j + 3): 



GQ v a (x,y) = ^- GO v a (x,y\e) le=0 = ~ 
ae ae 



( ,y» re* ^ 2-(-€ 

[ ^J—Q G v a (y,x\e) 



(yy) 



l+e 



+ 

|e=0 



= - < ^^ QG v a (y,x) + GQ v a (x,y)\n— + •••. (117) 

yy yy 

From this equation one can easily deduce that A G ®f a (x,y) = 0. The a and 6-kernels in the GG 
channel may be obtained form the QQ ones by means of the differential equation 

^- GG v i ( x , y ) = -2^%\x J y), for i = a, 6, (118) 
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which immediately provides the generic form (|62|) as a solution with A GG f l (x,y) = A GG f % {y). 
Since the lowest moment, i.e. J dx GG v l {x, y), is a constant, A GG f l (y) with i = a,b itself has to 
be a constant, which turns out to be zero. 

Now we come to the dotted c-kernels. Since the eigenvalues of the LO kernel in the GQ channel 
are given by the a-kernels in the AA channels, i.e. GO v c . = ^% c -/2) = Vj/6 [see Eq. (|83D], we can 
derive the following two differential equations 

^ GQ v c (x,y) = GG v%x,y)+ GG v a (x,y), ^ GQ v c (x,y) = -2 ^v a (x, y) - QQ v a (x, y), (119) 

where the added diagonal terms on the r.h.s are fixed in the way that the equations are solvable. 
Taking the generic form (pDl), we obtain two differential equations for the addendum: 



-^A GQ f c (x,y) = - — (2x-3), -^A GQ f c (x,y)=x(2x-3)-6xx\n-. (120) 
oy y ox y 

The solution A G ®f c = x 2 (2x — 3)ln| is fixed up to a constant, which by comparison with the 
lowest conformal moments turns out to be zero. 

Since the c-kernels in the mixed channels are related by supersymmetry, we may use Eq. ( |117| ) 
do find the dotted c-kernel in the QG channel from the GQ one. The GG entry follows then by 
integration from 

QG v c (x,y) = ^ GG v c (x,y). (121) 

Since A GG f c (x,y) = 2^(y — x) has been obtained by a direct calculation of a graph with a 
quark bubble insertion [see conformal predictions (|63"D 1 in Ref. f3lfl , we can use Eq. ( |121| ) to get 
AQ G f c = -^(4x-5y + 2xy). 

Now we are ready to present all dotted kernels. For odd parity we introduce the dotted matrix 

VUx ) = (° F [ QQ6 ^] + 2T FN f QG v(x,y)\ ^ I=u . (U2) 
{X,V> \C F G %\x,y) C A GG v\x,y) ) ' 

Here we use the decomposition as for the LO kernels ([74]) and (|75"D, 

QQ v = QQ v a + QO v b , QG v A = - QG v a , GQ v A = GQ v a , GG v A = \2 GG v a + GG v b ] , (123) 

QG- V _ QG- A _ 2 QG- c> GB^V _ GB-A + for g = {Q ^ Q} ? (m) 

and included the same "-(-"-prescription, although the dotted kernels are regular at the point 
x = y. The general structure of AB v % are given in Eq. (|60D, where the addenda are 

A AB f a = A AA f b = 0, A GG r{x,y) = 2- 2 {y-x) (125) 

A QGfc = / 4x _ 5y + 2 s A GQjc = x 2( 2x _ 3) ^ £_ 

3y y 
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3.2 Construction of G kernels. 



The construction of the diagonal G(x,y) kernel, related to the crossed ladder diagrams, goes in 
the same manner as demonstrated in subsection [2.3| for the reconstruction of the LO kernels from 
the QQ one by means of conformal and supersymmetric constraints. Since not all steps are quite 
obvious, we give here the construction in detail. 

The colour structure of the entries in the G kernel is 



G\x,y) 



l( 2C F (C F - c f)[^G I (x,y)} 2C A T F N f Q G G I (x,y) 



C F C A G QG\x,y) 



GG G (x, y) 



(126) 



for / = {A, V}. As explained in subsection |2.4| , from the result (|9~I|) - (j9~3|) in the flavor non-singlet 
sector [15|, [16], [H]] and its correspondence to the forward case given in Eqs. ([?|)-([56|), we conclude 
that all entries in this matrix have the generic form 



AB I 



S G' (x, y) = 6(y - x) ( A V + A^iT j y) + 6(y - x) ( A *W + A A *W) (x, y), (127) 
with the following expressions for H and H 

(128) 



(ABTj 1 



AB H\x,y) = 2 



± AB / (Li 2 (x) + In y InS) - AB f Li 2 (y) 



AB H\x,y) = 2 



Lio 1 



x 

y, 



+ -Wy\+ AB f {U 2 {y)-U 2 {x)-hiyhix) 



(129) 



Here the upper (lower) sign corresponds to the A = B (A ^ B) channels. In the QQ sector we 
have /S9®H = A^^ff = 0. However, in general it turns out that the non- vanishing addenda are 
needed to ensure the diagonal form of the kernels. The forward limit of these AB G ! kernels has as 
a generalization of Eqs. (|94T)-(|96D the following form: 



ab g\z) = \LlM AB G\x,y) 



z<0=>z>0 



B{z)9{l-z) 



AB H I (z)± AB H I (-z) 



for 



I = A 

J = V ' 

(130) 

where we employed the (anti) symmetry of the singlet parton densities to map the mixing between 



partons and anti-partons, given for — 1 < z < into the region < z < 1. The AB H f and AB H I 
functions are defined as 

AB H\z) = LIM( A V + A AB H I )(x,y)= AB p\z) (in 2 z - 2((2)) + A V \z), (131) 
AB E\z) ee LIM^ 3 !! 1 + A AB H I ) (x,y) = 2 AB V \z)S 2 (-z) + AB T\-z), (132) 



where we have to require that AB T and AB T are rational functions and/or terms containing single 
logs of momentum fractions. Note that the mapping, which took place in Eq. ( |130| ), can be also 
performed in the ER-BL representation by the substitution y — > y and is discussed in section f|. 



28 



First we reconstruct the missing three entries of A AB H A and A AE H A in the axial- vector case. 
Since the AB G A are essentially determined by the LO functions, it is convenient to decompose 
them in the same way as the LO kernels in Eqs. ( [72] ) and (|73|): 

QQqA _ QQqo, _j_ QQQb QGqA _ _QGqcl GQqA _ GQQa gg Q a = 2 GG Q a _|_ GG Q b (133) 

As for the kernels AB v A (x, y) to LO, we can obtain the QG entry from the a-kernel in the QQ 
channel by a derivative w.r.t. y as in Eq. flTop . Since the Q^G" kernel is defined in Eqs. ( |127| )- 
( p9|) with A QQ # a = A QQ ZT = 0, we find the addenda in the QG channel from 

A^ G H A (x,y) = - QG H A (x,y) + ^H a (x,y), A^&y) = - QG H A (x,y) + ^- QQ H\x,y). 

ay ay 

Since conformal covariance and supersymmetry connect these functions in a simple way with the 
GQ ones [see Eq. (|8|)], 



OQnA, , ^? Q GnA ( ^ ( A^H A (x,y) = -^ G H A (y,x) 

GQ G A (x,y) = - — - QG G A (y,x) { , ._ * v A 

W {AC ( 3H A (x,y) = ^Q G H A (y,x) 

we can write our findings in the following symmetric manner 

A^ G H A (x,y) = A^ G H A (x,y), A GG H A (x, y) = ^A GG JI A (y, x), (134) 
A GQ H A (x,y) = -A GQ H A (x, y), A GQ H A (x, y) = -2— \nx + 2— \ny. (135) 

y y 

To find the remaining GG entries we employ the differential equation (^), i.e. 





, QQ H A (x', y) + Z QG H A (x', y) + 3A® G H A (x', y) 
dy 



A GG H A (x,y) = - GG H A (x,y) - J* da/ 

and an analogous equation for A GG H A (x,y). Again the remaining freedom can be fixed by the 
necessary condition (|30D for diagonality w.r.t. Gegenbauer polynomials and from the requirement 
that the moments GG G A 1 vanishes for j > 1. To simplify the result, we remove a symmetric 
function^ (w.r.t. the simultaneous interchange x — ► x and y — > y) which enters in both A GG H A 
and A GG H A kernels, however, with different overall signs and, therefore, disappears from GG G A : 

A GG H A (x,y) = -A GG H A (x,y), (136) 

aGgtj a ( \ X ^ ~ x (x + y-3xy) x(x + y-3xy) 

A H (x,y) = - 2— - 2 — ma: -2 — \ny. 

y yy y y vv 

Now we come to the parity even sector. Instead of dealing with the whole sector, we can 
consider only the difference between vector and axial-vector kernels 

AB G V = AB G A + AB G S ^ Q G( f = _ 2 Q G Q^ GBtf = tfBg for B = { Q , G } , (137) 



J Here and below for the vector case this functions slightly differ from that one in Rcf. 
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and the same notation for AB H, AB H and their addenda. 

In LO we were able to write down the set (|7ID of differential equations that determine the 
entry in the GQ channel depending on the a-kernels. Unfortunately, the analogous equations for 
AB G 5 provide us a solution, which does not preserve the generic form of GO: G 5 in the forward limit 
given in Eq. ( |130| ) - ( |132| ). However, it turns out that we can restore the generic form by adding 
convolutions of c- kernels of the diagonal channels given in Eqs. ([74]) and (|77|) : 

f) r i 

— GQ G 5 (x,y) = -4 QQ G a (x,y) + 9 QQ v c ® QQ v c (x,y) - 9 QQ v c ® QO v c (x,y) , (138) 

OX I J 

^- GQ G s (x,y) = 2^ GG G a + 2 GG v c ® GG v c (x,y) + 2 G % c ® GG v c (x,y)y (139) 

where the kernel GG G a are part of the whole parity odd functions derived in the fashion already 
explained above. Note that this set of differential equations represents after separation of the 
^-structures in fact two sets, one for /SP G H (x, y) and the other one for A QG H (x,y). The two 
integration constants can be easily determined from the vanishing of the conformal moments 
G ®G c jo = for j > 0. Finally, we simplify the solution by adding pure diagonal pieces containing 
a and c kernels and their convolution as well as by removing symmetric terms that die out in GO: G. 
Using again the supersymmetric relation (|68| ), we write our findings in the mixed channels in the 
following way: 

A QG H\x,y) = -^A^lp), A^ G H 5 (x,y) = 7 %A G m S (y,x), (140) 



(yy) 2 ....... ( ?/?y ) 



A G V(x,y) = A G( m S (x, y) + 20 ^ ^ - A ^ 3 + ^ Ins + 4 X(3 + 2%) lny 

3y 3y 3y 

x 61 / 25 \ 

A GQ H (x,y) = 5 + 2x5(1 + — (a: -5) -(3 -105) In y + (3 -10a:) Ins) (141) 

9 V 18 / 

x (6 — 195 + 6a; 2 ) n x (y + x(x — x)) , x (y + x(x — x)) , 

H — - 2— = — hiy + 2— — lnx, 

3y y y 

The last task is to construct the GG kernel by means of the constraint ( |69"D where Q®G C = 0. As 

already described, the solution is obtained in a straightforward manner and reads: 

Kan^s, \ *artTT$, s x(3 — 135) 2a; (2 + 3a;) „xf„x — x 2 + 35\, 
A GG H (x,y) = A GG H(x,y) + — '- — ^-2- 2^— + ln5 

v yy v \ v y J 

x ( x — 5 2 + 3a; \ , , , 

-2- 2 + — lny, (142) 



v v y y 

rrTTS/ (1 -x 2 )(l -215) 5(19(1 + x) -36a; 2 ) 25 3 

GG H (x, y) = ^+ — i H 

1 y) 3y 2 3y 3y 



( x 3 x 2 (21-20x) xx 2 \ fx 3 5 2 (21 - 205) xx 2 \ 

+2 — i - - 2— )\nx + 2 — i = - - 2 lny. 

\3y 2 3y y J \3y 2 3y y J 

This completes our construction of the complete G matrices relevant for parity odd and even 
sector. 
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3.3 Restoration of remaining diagonal terms. 

Now we have to find the remaining diagonal pieces of the ER-BL kernels, which are expected to 
have a simple representation in terms of the kernels known from Eqs. (|7^)-(|7H|). The DGLAP rep- 
resentations of the missing terms are obtained with the help of Eq. ( |115| ), where we can interchange 
the exclusive convolution with the forward limit procedure: 



D\z) = P {l)I (z) 



-P^tp^ + ^t 



9 



k)P {0)I 



G 1 > (z). 



(143) 



Here we have introduced the analogous notation for the DGLAP representations as used in the 
non-forward kinematics. The DGLAP kernels to LO read 

for / = {A,V}, (144) 



P^ 1 = UMV^ 1 



F l\ f ^{z) 



C F GQ p\z) 



Y{z) -f5(l-z 



where the entries are decomposed in the same manner as in Eqs. ( |74"D and (|75|) , i.e. 



QQ _ QQ a , QQ b QG A _ _QG a GQ A _ GQ a GG A _ r, GG a , GG b 

P — P * P P — Pi P — Pi P — * P > P 

QG p V ^QG p A_ 2 QG p c^ GB p V ^GB p A + 2 GB p c for B = { Q^ h 



(145) 

(146) 



and the functions AB p l are defined in the following way 



1, 



2z 



QQ^fl — _GG„a -| ^ QQ^b _ GG^b _ _Z_ QQ,/- ± ( ' r V ' ' (147) 



-on" — _~~Vi a — 1 — r _ ^^rP 

P — o P ~ 1 ^1 P ~ P -| 

2 1 — z 



P = 3 V 



:i-^) 3 

3z ' 



QG a 



■p a = 1 - 2z, QG p c 



:i-z) 



GQ a 



<p a = 2- z, 



GQ p c 



:i-z) 



The "-(-"-prescription is now uniquely defined in the QQ and GG channel as the conventional one: 



Mzj\, = M p(z)-8(l-z) / dy AA p(y). 



The limit of the dotted kernels defined in Eqs. ( |122| )- (|125|) give 



P 1 ee LIMV 1 



C F Q®p(z) 2T F N f ^(z) 
C F GC tf(z) Ca^p^z) 



(148) 



with 



Ap _ AApA ln z _|_ A^a 



^V = AB p A lnzT QQ p a for 



AS = 

AB = GQ ' 



Q G p v = QG p v(\nz + ^)- ^ GQ p c + l ^ QG p a - Q V, GG: p v = - G V, 



GG p V _ GQaA , GG„c 



11 



GG^c 



'p A + 2 Lr V In z H ~~jr. 

3 



(149) 
(150) 
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The limit of the g matrix defined in Eq. ([58]) is 

I -C F 



g(z)=UMg(x,y) = 2 



ln(l-z) 
1-z 



V 



ln(l-z) 
z(l-z) 



) 



(151) 



Finally, we need the complete forward limit of the G kernel given in Eqs. ( |130| ) - ( |132j ), i.e. we 
have to determine AB T I and AW T 1 . From the result of the G matrix given in Eqs. (|126|) - (|129|) 
with the addenda dl3l-(p6|) and (|T40|)- (p6|) , we obtain the following expressions for the sum or 
differences in the parity odd and even sector, respectively: 



AA T A {z) + AA T A (-z) 



4 A Y(-z) ]nz + 8 A y< 



iB T A {z) + AB T A (-z) = T^ BA p a (-z) In z T l2 QQ p a {z) for 



AB = QG 
AB = GQ 



QQr v (z)- QQ T v (-z) = o, GO r v (z) 



G 



■- GG p c (-z)\nz, 



(152) 



(153) 



iB r v ( z ) - BA r v (-z) 



8 -^(~z)(2±3lnz) + f^(z)\nz for l^Z% 



With these results we are now able to find D I {z) from the known DGLAP kernels in NLO by 
means of Eq. (|143|) . It remains a simple exercise to express these findings in terms of convolutions 
of a,b and c-type kernels defined in Eqs. (|147 ). In the following two subsubsections we treat the 
parity odd and even sector separately. We start with the parity odd sector responsible for the 
evolution of axial-vector distribution amplitudes, since the structure is simpler than the one in 
the parity even sector. 



3.3.1 Parity odd sector. 



We already have explicit expressions for all kernels entering on the r.h.s. of Eq. ( 143 ). It remains 
to convolute the P and g kernels with the LO ones. Since cancellations between these separate 
terms drastically simplify the final result, we present it as entries of the matrix 

ji _ p 1 X I n (nu , Pa 



p(0)I + M> j 



9 



(154) 



In other words this matrix contains the diagonal part that we could not separate from the confor- 
mal anomalies. The calculation is straightforward and results in 



qq a a = 2C F T F N f {4(1 - z) - (1 - 3z) In z - (1 + z) In 2 z] + N! U, 



(155) 



where the convolution in the non-singlet part NS ^4 = —C\ ^^p^^p — Cf^q^^p/2 has been worked 
out in Eq. (|107|) . The remaining channels read: 

7 1 



QG A J 



2C F T F N f \ (l-2z) 



-- + 2C(2) - -\n z z + 2\nz\n(l - z) - z) 
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+ (1 - z) [-4 - 3 lnz + 41n(l - z)\ 



+2C A T F N f \ (1 - 2z) 2C(2) + ln 2 (l - z) - 2(1 - z) [1 + 2 ln(l - z)] 



(156) 



GQ^ 



-61nz - 2(1 + z) \n 2 z 



d \ + In 2 - 41n(l - z) + (2 - z) 



ln 2 z - ln 2 (l - z^ 



A) 



-C F y {3 - ^ + (2 - z) In z] + C F CV - 5(1 - z) + (4 - 9z) lnz 



(157) 



3 - 2z + 3z 2 



ln(l - z) + (4 + z) In 2 z + (2 - z) ln(l - z) [ln(l - z) - 2 In z] 



GG J.A 



A 1 



; 2-3z + 2z 2 
1 - z 



{-2C(2) + In 2 z - 2 In z ln(l - z)} 



(158) 



+4(1 - z) [4 - ln(l - z)] + 12(2 - z) In z + 4(1 + z) In 2 z 

- + 

+2C pTpNf {-10(1 - z) - (7 + z) lnz - (1 + z) In 2 z} . 

Now we are ready to extract the remaining diagonal terms AB D A from 

D A (z) = P^ A (z)-{A A + G A }(z). (159) 

As observed in the flavor non-singlet case, all double logs and Spence functions appearing in the 
DGLAP kernel in NLO g§ |47| are contained in the sum of A A given in Eqs. (p^55|)-([T58|) and G A 
defined in Eqs. (poD -fpg) together with Eq. ([TBI) . The remaining single log terms give us a hint 
to write the entries of D A (z) as convolution of the simple kernels (|147| ) , which allows us to restore 
the ER-BL representation. In the following we present this issue in detail. 

In the QQ channel the only new information arises from the pure singlet term. Employing 
formula (|159|), we immediately find that this term is 



-QC F T F N f QQ p a -QC F T F N f QQ u a . 
Adding the result of the flavor non-singlet sector, the QQ entry reads 
QQ D A = ns d _QC F T F N f QQ v a , 
NS D = Cl\D F ], -C F ^\D ; ] -C,[C 



c F % [DA 



+ 2 QQ v a 



(160) 
(161) 



where D F , Dp are given by Eqs. (|112|) and ( |102| ), respectively. 
In the QG channel we find in an analogous way 



QG D A (z) 



3 C F T F N 



l + 2z) lnz - 3(1 -z) 



:i-2z, 



-2 C A T F N f {-3(1 + 2z) In z - 9(1 - z) + [1 + 2((2)} (1 - 2z)} 



33 



by a convolution of a-kernels, 

QQ p a^QG p a = _Q + 2 z) In Z - 3(1 - 

together with ® G p a = 1 — 2z the ER-BL-representation: 

QG D A = 3C F T F N f (^v a ^ G v a -^ G v a \ 

-2 C A T F N f (3 QQ v a ® QG v a + [1 + 2C(2)] QG v a \ . 



(162) 



A more complicated expression arises in the GQ channel: 

GQ D A (z) = (7f,{(2-z)ln(l-z)-(4 + z)lnz-6(l-z)-^(2-z) 

-C F ^( 2 (2-z) ln(l-z)-(2-z) h, z + %z-\{2-z\ 
2 1 2d 



-C f Ca 



( 3 + 2z - z 2 



V 



ln(l - z) - 21nz - 6(1 - z) 



2C(2) 



(2-,) 



which is characterized by new terms containing (2— z) ln(l—z) and [ln(l— z)]/z. Such contributions 
can be generated by convolutions of 6-kernels with a- and c-kernels, respectively^ 



G V ® r 



P 



2(2 - 2) ln(l - z) - (2 - z) In z + 3/2z, 

GG p A ® G V = 2(2 - 2) ln(l - z) - 2(4 + z) In z - 12(1 - z), 
(1 - z) 2 

ln(l - z) + 2(4 + z) In z + 10(1 - z) 



G V® G V = 2 



GG p a ® GQ P a -2(2 + z)lnz-6(l-z), GQ p a = 2- 



z. 



The result can be written as 

1 



gq d j 



(% \ ^ GG v A ® G( V - ^ G Vl - ^ G( ->r" 



_c F c A r- GG v A ® G( ^v c + 



2 GG i) A - - GG v 



h 
2 



» G V- 



GQ 



3 " 2C(2) 



The remaining entry reads in the DGLAP representation 



'D (z) = C 2 |4(l-z)ln(l-z) -2(5 + 3z)lnz + 

-25(1 - z)} - |2(1 + z) lnz + H 

— C F T F Nf {-4(1 + z) In 2 - 10(1 - z) + 5(1 - z)} 



2 - 3z + 2z 2 
1 - z 
2 - 3z + 2z 2 
1-z 



(163) 



2 v 1 



6(1 - z) + 28(1 - z) 



3 Note that we can eliminate one of them e.g. GQ p a ® [ QQ p] = [ GG p A ] ® G Qp a - 3 GG p a <3 GQ p a /2 + 3 GQ p a /2. 



34 



and can be expressed by convolutions of a- and 6-type kernels, 

GG p A^GG p a = ^ _ ^ ^ _ ^ _ 4(2 + z) hi Z - 16(1 - z) 

2 — 3z 4 

'p~ = 2 ' , ^ u p a = 2(1 - z) 

from which the desired kernel follows: 



G V® G V = -4(l + z)ln Z -8(l-4 G V = 2 2 * Z _ +2Z2 ' GG ~ a 



GG D A 



c 2 



GG A l GGv a 

2 



r,G v a 4 - GG v A - - GG v a - 25(x - y) 



2 

3 " I 

00 f 1 GG„. a X GG„. a i ^GG„A GG„, a 



- C F T F N f \ GG v a ® GG v a - GG v a + 6(x-y) 



(164) 



These results provide us the missing information to obtain the whole kernel from Eq. (|114 ) in the 
flavor singlet sector for odd parity. 



3.3.2 Parity even sector. 

In analogous way we treat now the parity even sector, starting with the calculation of the matrix 
«4. v . For brevity we present here only the results for the difference between vector and axial- 
vector case, i.e. A 6 = — A. defined in Eq. ( J154 ) . This difference is induced by the c-kernels. 
Thus, it is not surprising that ® G A S consists of terms \ jz In* z In 7 ( 1 — z) and z 2 In* z In- 7 ( 1 — z) with 
< i + j < 2: 



QQj8 
QGj& 



GQaS 



-2C F T F N f |— (1 - z)(31 4 142z - 5z 2 
2C pTpNA (1 - z) (10 - 7z) 4 2(1 - zf 



3z 



'l + 7z+ 7z 2 - z 3 ) Inz 



-2C(2) + 31n- 4 In 



1 - z 



1-z 



-(1 - 4z 4 (1 - 2z) In z) mz j + 2C A T F N f ^ 
— -(4 4 13z4 88z 2 ) Inz 4 (1 - 2z) In 2 z 
-(1- z? Uc,{2) 4 \ ln(l -z) + 2 ln 2 (l - z) 



1-z 
18z ' 



C 2 F \ -3(1 -z)- [32 4 (2 - z) Inz] Inz - 2 



'1-z) 



z 



ln 2 (l -z) 



+C F C A \ -(17 4 3z - 21z 2 4 z 3 ) + —(31 4 24z 4 57z 2 - 4z 3 ) Inz 

1 2(1 - z) 2 

4-(2 - 2z 4 z 2 ) In 2 z 4 V - ; [-3 - 2 Inz 4 ln(l - z)\ ln(l - z) 

4C F /3 {1 - z 4 (2 - z) Inz}, 



(165) 



88 4 175z 4 367z 2 ) (166) 



(167) 



35 



GG aS 



C 2 



2<yl Z \ ll2 - 2z + l\2z 2 ) + ^-(22 - 18z + 81^ 2 - II2 3 ) lnz 
9z 3z 



1 y , 

2(1 -2) 3 



3z 



6C(2) + 11 ln(l - z) - 3 In 2 2 + 6 In z ln(l - z) 



(168) 



+2C F T F N f {4(2 - z)(l - z) + 2(3 - z) \nz} 



Taking these results and the corresponding ones for the axial- vector case in Eqs. ( |155| )-( |l58| ) 
as well as the definitions ( |130| )-( |T3~2D together with ( |153| ) for the G kernels, we obtain from the 
NLO kernel EH] the diagonal terms 



D v (z) 



P^ v (z)-{A A + A 5 + G v }(z). 



(169) 



Again we observe a cancellation of all double log terms. As to be expected, the convolution 
formulae from the previous subsection are not able to express all terms in D v as convolutions. 
What is missing are convolutions of c-type kernels with a, b and c ones. We show in the following, 
that this is enough to restore D v for non-forward kinematics. 
For the pure singlet part in the QQ channel we have 



2C F T F N f J - QQ p a ® Q V - 2 Q V + 



17 + A6z + 17z 2 AM .1 + 8^ + z 2 , 1 
h4(l + z) lnz 



3z 3z 
where the convolution of a kernels is QQp a ® QQp a = —2^p a {z) — ®Qp a (—z) In 2. Making use of 

11 + 38z + llz 2 



QQ p c g, QQ pC 



1-2 



27z 



.I + 8Z + 2 2 OG . (I-2) 3 

I + 2) : Inz, QQ p = 



9z 



3z 



we can restore the ER-BL representation 



QQjy 



NS 



'D + ACpTpNfi - QQ v 



c ® QQ v° 



a 7 

+ 6 



(170) 



For the QG entry we find 



Q G D V (z) = 2C F T F N f \ - - QQ p a ® Q V - - V - 3 Q( -' 

I 2 4 



'1 - 2z + 2z 2 )lnz 



+2(1- z) 2 ln(l- z)\ + 



2C A T F N f { — QQ p a ® Q G p a + f ^ 



+ 4C(2)) Q V + 8(1 



301 



+^(l-z) 2 ln(l-z)|. 



\ 3 

7 + 130z + 70^ 2 
' 9~2 



9 



2C(2) 



QGpa 



+ -(20 + 44z + 5,2 2 )lnz 
3 



With the help of 



P 



k> QG p C 
Q p c ® QG p° 



1 - z + (1 - 2z + 2z 2 ) In 2; — 2(1 — z) ln(l — z), 

s 1 + 19s + IO2 2 1 , , 
-(1-*) ^ -(3 + 6^ + z 2 )ln^, 
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we can express the remaining terms as convolutions and obtain so the ER-BL representation: 



QGjf = -C F T F N f h 
+ 2C A T F N f i- 

+ 



, Q%° + QQy 0, d Q V + - QG v a + Q QG v 



+ 56 w v 



130, 



f-2C(2) 



QG v a 



301 



'v + 



4C(2) 



171^ 



A similar situation appears in the GQ and GG channel. Beside the known convolutions we need 
only the following ones 

fl - z 2 ) 

+ (4 + z) In z + ^ ^ ln(l - z) , 





® GQ p C 


= 5(1- 


GG / 


h GQ p c 


= -(1- 


G y 


h GG p c 


= 10(1- 


GG p c 


h GG p c 


= "(1- 



10 + 19^ + z 2 1 + 6z + 3z 2 
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11 + 38z + llz 2 



In z, 



(l-z) J 



hi(l-z), 



3z 



1 + 8Z + Z 2 
1 + z) mz, 



and thus find the following diagonal kernels 



GQjy 



gg A e 



-C F /3 { G S 
+C F C A \ 



A e 



- GG v a + 3 G( V 

[.GQy 11 _|_ GQy C 



- 5 GG v a <g> GQ y a - 3 G( V j 

+ - GG t> <g> Giy i> + - H; ^ 
4 3 J 



GG V (g) 



2 



—GG v a (g, G( 3y a _|_ 9 G( V (§> GQy c 
6 



/43 
V~9~ 



+ 2C(2) G V+ ^-4C(2) 



9 



GQ,- 



C 



2 J GG 



A e 



(g> 



LrMy _L y 

3 



14 G( V <g> GG t> a + 12 G( V (g) GG v° 



+ - GG v 



A 131 



GG„, 



91 



12 + -"V - 25(x - 



C 4 —— s — — V 



(172) 



(173) 



oa « A GG^ a + 5 GGu v + 3 g V + 13^ + 2J(x _ ^ 
3 3 

+ C F T F NA GG v a ® GG v a + ^ GG v c -5{x-y)Y 

It is worth mentioning that our result for the evolution kernels in the parity even singlet sector 
possesses the correct conformal moments in both the physical and unphysical sectors. This is 
to be contrasted with an explicit momentum fraction space calculation at LO and quark bubble 



insertions in NLO kernels for the mixed channels 2S where the improved kernels do not appear. 
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4 Explicit representations of the kernels. 



In this section we give the explicit form of the ER-BL type kernels. The analytical continuation 
to the whole (x, y)-plane allows to extend our results to the skewed kinematics and extract from 
them the so-called skewed DGLAP kernels, which are responsible for the evolution of SPDs in the 
DGLAP region. The exclusive convolutions we need are straightforward to handle and details are 
collected in appendix |C|. To represent the final result in the shortest possible manner, we expanded 
the output of the convolutions in terms of powers of logs. Then we rewrite their arguments in 
terms of the ratio x/y where the remaining ones posses only an x-dependence. It turns out that 
these terms for the (non-)diagonal entries are (anti) symmetric with respect to the interchanges 
x — > x, y — > y, so they have the support < x, y < 1. 



Now we should comment on the G-kernels constructed in section [3.2| . Although they are 
diagonal, this property is lost for the unphysical moments when we map the second ^-structure to 
the first one by means of symmetry. A slight change in the definition of AB H J and AB H I in Eqs. 
and ( |129| ), by subtracting a symmetric function, i.e. 



2 [6(y -x)± 6(x - y)\ 



AB fU 2 (x) ± s 

from the common ^-structure and adding them to the second one, ensures that both parts are 
separately diagonal. Thus, in the following we use instead of the definitions ( |128| ) and ( |129| ) the 
definitions: 



AB h J 
AB h l 



±2 AB f ! In x In y - 2 AU f [Li 2 (x) + Li 2 



ABrl 



;i74) 



ABrl 



2Li 2 



l- X - 

y, 



In y 



+ 2 AB f [Li 2 (y) - In x In y] ± 2 AB ~f Li 2 (x) , 



ABJ 1 ! 



where again the upper (lower) sign corresponds to the (non) diagonal entries. Together with the 
addenda given in eqs. ( |134|) - (|136| ) as well as eqs. ( |14q) - (|142"|) the kernels are separately diagonal 
w.r.t. the Gegenbauer polynomials. 

The LO kernels are defined in Eqs. ([7^)- ([78]) and the NLO kernels we write for convenience as 



V WI (x,y) 



+ 





o GG vS ]I 



AA 



v® 1 



where contrary to our previous definitions 
( |79|) applied also in the GG channel. The contributions concentrated in x 
the lowest conformal moment. We have in the GG channel: 

N 



5(x-y), (175) 
denotes now the conventional regularization 



y are simply fixed by 



GG,AV)I 

v ll 



GGy(A)I 



108 



(35C A + 74CV) 



(176) 



95 
27 



14 



N 



C(2)+2C(3))+^(29Cx-28C. 
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The QQ channel contains a contribution coming from the non-singlet kernel and the so-called 
pure singlet part. Analogous to the forward case we define the ± kernels in the flavor non-singlet 
sector to be 



cUe( y -x ! 1 



2C(2)) + - (^ QQ f - |) In X - - (QQf - In 



x {5QQ f+ * + QQ f l n *)_C F ( C 
13 y y) 

T QQ h(x,y)\ + 



x 

c A 



— lnx(l + lnx-21nx)| - C F ^-6(y - x) (177) 

%-*)(! W / + 2- + < *k(s,i/) 
[3 y 



a; — > a; 

y^y . 



+ 



1 ± 1 



C F Cf — 



13 



-6C(2)+4C(3) 



5(x - y). 



The kernels QQ V ( ' are responsible for the evolution of the sum (difference) and difference (sum) 
of quark and anti-quark distributions in the parity even (odd) sector. Our results for the entries 
of the parity odd sector read: 



QQy {1)A = 

QOyW- _ 2C F T F N f [e{y - x)- In - ( 1 - In 

V V 



- --lnx(l-lnx) + 

yj y [y-+y 



QGyi 1 



)A 



(178) 



(179) 



N f T F C F \e{y - x) [2 (5 - 2C(2)) QG f A + — +( 2^ G f A - % + —) In - - 2^ G f A + ® G J A 

I L yy \ y yyj y v 

+1) i n f 1 _ ln2 x + , QG A + Q ^Ax ln2 /y _ x 

yy J \ y y v y Va; 



- QG / A lnx(2 + 21nx 



3a; 



-lnx) + — ma^+AyrFCJ^y-a;) 



In 1 



a; 



yy, 



+ QGjA 



3a; 



—® G h A (x, y) — ( ^ G h, A (x, y) 



+ 



6xy — 4xy 

yy 2 



_ 4 QG f A _^_ 2 (QG f A_2A ln X +2 ( QGjj 

yy \ y J y \ 



3a; 



lnx + °^f 4 lna;(2-lna;) + — hvM - 



^fUy-x) 



x — > x I 

y^yj ' 

(180) 



3 ^(l + g) + a: 2 ?/ + 2x(x - a;)y ^ x _ x(x + y) + a;y ^ / _ x\ + gq^a 

y yy y yy \ 



y, 



X 



In 2 — — ^GQjA _|_ GQjA^j x _ 



X 



a; 2 y + (1 - a; 2 )y 



lnx- GQ / A ln 2 a; - C, 



A) 



x<^(y-a?) 



_ 2 GQf A + 2 £ + ( GQ/ A + GQjA^ ^ ^ _ ^ + GQj-A ^ + 
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x{9(y-x) 



-2a2)) GQ f A --- 



' x[x — x) xx 



X 



+ 5— I In - + 

y y J y 



x — 2x 2 x — 2x 2 
+ 



y 



y 



X 



b t 1 - ~\ ( GWjA + 3 j) '" 2 ? + 5 ( CQfA + G? ") '" 2 (i - - 



+i G % 4 (x, !/ ) + 



, JU JU JU 



y y 



GQ f A (-2 + -\nx-\n 



x 



3^ 
2? 



In a; In x 



x — > x \ 

y^y}' 



GC V (1)A = 

C 2 A U{y-x) 



(181) 



^ - 2C(2) ) GG f A - ^ - 6- 



15 x 2 p xx t x(x + y) 



yy 



y 2 y 



IX XX X 

- 2--6— + - 

V yi/ r y 2 



X 



In - + f GG f A + 2 X -\ In 2 * - ( GG f A - ^ ) In X - In f 1 - ^ - i GG //(x, y) 

y V y 1 ) y K ' y \ yj 2 



+C^(y - x) 



( (x — x)x xx x x , , 3x — y . _\ , 
+ 3 = o— + — — 2— (Inx + In an H — — mx ma: 

V 1/1/ v v y y 



3 J 2 y 2 



5x 2 x(y + xy) x 



x 2 



y 2 y 



+ — lny + — lnx } + C F N f T F \ 6{y - x) 



x 



x(l — 2xy) x 2 , x x , x / , a; 
2 — — - — — + 2— In - - — In- 2 + xln- 

yy y y y \ 



y 2 y 



v. 



X 



+ — In x (2 + x In x) > + 

y J 



a; — > x 

y^y 



Except of the pure singlet part, where the convolution of two c-kernels generates Spence functions 
which die out in the forward limit, we recover in the parity even sector a similar structure as 
observed in the previous cases: 



QQyW = 
QQyW+ + 2C F T F N f \e(y - x) 



(182) 



x(3 — 8xy) x(5 — 8x) , x fx 
— + In - + - - 4xx 

y y y \y j 



290 



x In 2 — h 8xx (Li 2 (x) — Li 2 (y) + lnx lny) — x v 

y ; 

x In x — f — — 4xx ] In 2 x + 4xx In x In x \ + < 

\y / J I?/ — ► 3/. 



'x(5-8x) 2x(9-19x)' 
, V + 3 . 



QGyW = 

N f T F C F \e(y-x) 



(183) 



2(5-2C(2)) 



x 



ln X_ 2 (QG f V + QGjV + ± 



yy, 



; f v — ^QGj-A _|_ _|_ 2 I QGyv _|_ QGJ V _|_ 2^ 2y + xy 
yy V 2 ?/y 2 
+ ln2 x + QGfV + QG ^v ln2 /y _ A 
y Vx / 



In 1 - 



-2 (® G f v + ^ + ^ + X( ^_y 2) ) ^ x - 2« G / y In x In x + ( QG / y - QG f V - QG f A ) 



x In 2 x j + jV/ZVCUjflfa - x) 



6x(3-4x) 2x(ll-16x) 4x(l - 3x) 



+ 
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Now we will derive the NLO skewed DGLAP kernels which govern the evolution of SPDs in 
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the kinematic range of ( < z < 1 and — 1 + ( < z < with z = and ( = We chose this 
set of variables such that we have the closest resemblance to the usual DGLAP kinematics. In 
fact for DVCS and vector meson production ( = x B j- 

After the continuation of the ER-BL kernels to the whole x — y plane by the replacement of 
the ^-function structure, we derive the kernels in the following way which is suggested by the 
non-zero support of the ^-functions after the replacement x — > z/(, x 
y — > 1 — l/£ = ^ for the above mentioned kinematical regime: 

QQfWi (1^ _!_ (QGf(i)i + QGj(i)J 1 1 I 



l-z/(Ef,y^l/( and 



P (1)I (zX) 



l (QQf{l)I ^ l 



GQj(l)/ ^ 1^ + GQj(l)J Q ? 1^ 1 (GGj(l)I Q ? 1^ _ GGy(l)/ 1^ 



, (186) 



where the functions AB f^ I (x,y) are the coefficients of the 9(y — x)-term appearing in the ER- 
BL kernels. Note that we omitted the additional factors of 1/z in the GG and GQ kernels so 
as to take into account the fact that in the forward limit, the skewed gluon distribution turns 
into zG(z,Q 2 ) rather than G(z,Q 2 ). Thus, we can take the "-(-"-prescription from the ER-BL 
kernels, where contributions concentrated in x = y turns into end-point concentrated ones, i.e. 
5(x — y) — > 5(1 — z). 

Following our afore mentioned prescription, the result for the parity odd and even skewed 
DGLAP kernels are up to end-point concentrated terms, which can be easily restored, given by: 
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where the and are given by 
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The LO skewed kernels appearing in the above formulas are given by: 
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where in comparison to our previous conventions we included here a factor 2 in the LO kernels of 
the QG-channel. The AB k\ were found to be: 
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The set of explicit formulae given in this section represents the main result of the present paper. 



5 Conclusions. 

In this paper we have given a detailed description of the formalism for construction of the two-loop 
exclusive evolution kernels in the flavor singlet sectors and presented explicit results. Our formal- 
ism allowed us to avoid complicated NLO calculations and was based on three main ingredients: 
(i) the known form of conformal symmetry breaking counterterms for renormalization of conformal 
operators in NLO transformed into the language of ER-BL kernels; (ii) supersymmetric constraints 
which allowed us to construct the contribution of cross ladder diagrams; (hi) reduction formulae 
and known two-loop splitting functions which completely constrained the diagonal part of NLO 
ER-BL kernels. The kernels predicted here are numerically checked by direct comparison of their 
Gegenbauer moments to the anomalous dimensions of conformal operators whose correctness was 
shown by supersymmetric |50| and super conformal |44] constraints. Moreover, the predicted /3q 



terms have been checked by a direct calculation of diagrams with quark bubble insertions. Note 
also that the predictions of the conformal operator product expansion rotated to the MS scheme 
coincide with the NLO coefficient functions. This proofs again our results of conformal symmetry 
breaking in QQ and QG channels. 

As a byproduct, our understanding of the general structure of ER-BL type kernels implies 
the discovery of a simple structure for the DGLAP kernels, which results from the topology 
of Feynman graphs. Indeed, the only new functions appearing in NLO is AB G(z) arising from 
the crossed ladder diagrams. Everything else can be represented as convolutions of LO kernels 
and the ones obtained in the forward limit of the conformal anomalies. This is interesting, but 
unfortunately not very restrictive, since the conformal algebra does not constrain these kernels. 

With the ER-BL kernels we have calculated, we open up a new way to effectively perform the 
evolution of singlet distribution amplitudes and skewed parton distributions by direct numerical 



47 



integration of two-loop evolution equations. The numerical algorithms are still to be developed 
but a priori it is clear that they will provide a superior alternative to the ones which rely on the 
orthogonal reconstruction approach. 

A.B. was supported by the Alexander von Humboldt foundation at the initial stages of the 
work. A.F. was supported by the E.U. contract # FMRX-CT98-0194. 



A Uniqueness of extension. 



In Eq. ( p3|) we gave a simple recipe for the extension of the ER-BL kernels into the whole region. 
Based on the holomorphic properties of the Fourier transform of "y(t, t', rj) we give here a more 
complicated method which proves the uniqueness of the procedure []J . Due to the scaling relation 
(|15D we can restrict ourselves to rj — 1. 

First we perform a Fourier transform of the anomalous dimension "y(i, i', 1) with respect to t, 
restricted to the ER-BL region \t,t'\ < 1. Due to the known support properties (see fig. |I]b) it is 
sufficient to restrict only the variable t', i.e. \t'\ < 1 (the anomalous dimension matrix vanishes 
then for \t\ > 1). Taking into account the support restrictions arising from the connection fll4] ) 
with the anomalous dimension of light-ray operators, we observe that 

7(VV|<i = |rfte iA ' 7 (t,0|i'|<i (199) 



«y(y, z) i\ QG i(y, z) \ iX{t , {1 _ y _ z)+y _ z} 



1 rl I Q 

o dy J dz \ 1.<£ZU "acj^y I 6 ' 



is an entire function in A and t' for \t'\ < 1. Thus we can perform the analytical continuation 
which coincide with the Fourier transform of "y(t, t', 1). This proves the uniqueness of the extension 
procedure. 



B Forward limit. 

We give now a more heuristic derivation of the forward limit procedure ( |20|) of the generalized 
ER-BL kernels. The moments of the kernels are related to the anomalous dimensions of local 
twist-2 operators by 

On the other hand the diagonal entries of the anomalous dimensions are given by the moments of 
the DGLAP kernels 

/ QQp QGp\ i /QQy QG \ 

dzz \ GQp ao P )^ = -- 2 {o2 i Gc 1 ) j - ( 201 ) 
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To extract the diagonal ones from Eq. (|200| ), we substitute y by 1/rj and multiply each entry 
on both sides with a sufficient power of 77. After extension of the kernels, we can rescale the 
integration variable and find in the limit rj — ► 0: 



■n^oj \-q\ \nGQy Iccy ) 1 77' 77/ 2 \ G Qj GG y 



(202) 



Interchanging integration and limit and comparison with Eq. ( |201| ) provides us with the desired 
formula: 

I ( QQy IQCV\ ext (z l\ 
P W = LIM^,,).lim-^ GQi/ | G(v j (-,-). (203) 

Let us look more closely at the forward limit of the +-prescription. In the QQ-channel the 
formal prescriptions for exclusive and inclusive channels are simple related to each other: 
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QQ V{x, yj\ = UM QQ V(x, y) - 5(1 - z) / dzLlM QQ V{x, y) = \ QQ P( 



(204) 



However, in the GG-channel an extra z factor appears which induce a finite term concentrated at 
z = 1: 
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V(x,y)-5(x-y) / dz^V(z,y) 
Jo 



z GG P(z) - 5(1 -z) J dzz GG P(z) (205) 
GG P(z)} + 5(1 - z) [ dz(l - z) GG P(z). 



In order to have the same simple correspondence between the +-prescriptions for exclusive and 
inclusive kernels as given in Eq. (|204|) for the QQ channel, we redefine, as implicitly done in our 



previous papers [p2 |, the standard +-definition by a finite part: 



GG] 



V(x,y) 



GG] 



V(x,y) -5(x-y) 



(206) 



dz GG V(z,y) + I dz(l - z) GG P(z) 
If the kernel v(x,y) contains only the usual ^-structure, i.e. 9(y — x), the limit yields 

UMv(x,y) = 6(1 -z)G(z)UMF(x,y), (207) 
UMF(x,y) = 



I / (QQp _ QQp^j l(QGp ^QGp^x / 

\ U(GQp + GQp^ l(GGp _ GGp^ ) ^ ) ' 



In the case of 6(y — x) -structure, we have to replace in the above equation 9(1 — z)6(z) by 
6(l + z)6(-z). 



C Exclusive convolutions. 

In this appendix we present the exclusive convolution formula. Due to the two different regions 
y > x and y < x appearing in the kernels, the convolution looks more complicated than for the 
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inclusive case. Fortunately, it is sufficient to consider one of the regions and for the convolution 
of two regular kernels, 

' x — > x ^ 



1)"W+"(B) ; - ' ~ K with i = {1,2}, 



AB v i (x,y) = 6(y-x) AB f(x,y)- 
we use their symmetry and define the convolution of the fi(x,y) functions by 
AB v 1 & BC v 2 (x,y)=0(y-x) AB f 1 & BC f 2 (x,y)-(-l) 



y^y 



(208) 
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_ ( _ ir (A) + ,(B) f* dz AB f ( 22) BC f ( y (21Q) 
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In the case that the kernels need a regularization, e.g. the conventional one 

(211) 



A "f4*,v)= A "f&v)-S(*-v) / dz AB f(z,y), 

Jo 

we have to arrange the integrals in such a way that each of them is defined. Since the regularization 
only appears in the diagonal channels, we have to treat the following two cases for convolutions 
in the mixed channels: 

AB fi® BB f2, + = [ y dz{f l (x,z)-f l (x,y)}f 2 (z,y)+ f dz[h{x,z) - f x {x,y)\ f 2 (z,y) 

J x Jy 

- [ X dz{f 1 (x,z) + f 1 (x,y)}f 2 (z,y), (212) 
Jo 

f AA ®f AB = f y dz[f 1 (x,z){f 2 (z,y)-f 2 (x,y)} + {f 1 (x,z)-f 1 (z,x)}f 2 (x,y)} 

J x 

f 1 dz[f 1 (x,z)f 2 (z,y) + f 1 (z,x)f 2 (x,y)} (213) 
J y 

+ / dz[f 1 (x,z){f 2 (z,y) - f2(x,y)} + {fi(x,z) - fi(z,x)} f 2 (x,y)], 
Jo 

where A ^ B and = {f AA ,f AB } and each term is separately integrable. The corresponding 
equations for the convolution of the diagonal kernels follows by changing the sign in the third or 
second integral, respectively. For the convolution of two kernels with "+" -prescription we write 

x — > x 



Vl 



v 2 



d{y-x) AA f l ® AA f 2 {x,y) + 



AAj 



y^y 



(214) 



+ 



where the convolution on the r.h.s. is given by 

ft A ® f2 A = I' dz{[f\(x,z) - h{x,y)} [f 2 (z,y) - f 2 (x,y)\ + [f\(x,z) - f\{z,x)} f 2 (x,y)} 

J X 

+ f 1 dz{[f 1 (x,z)-f 1 (x,y)]f 2 (z,y) - [h{z, x) - f ± (x, y)] f 2 (x, y)} (215) 

+ dz{ [A (x, z) - /i (x, y)\ {f 2 (z,y) - f 2 (x,y)] + [f l (x, z) - j x (z, x)] f 2 (x, y) } . 
J o 



50 



References 

[1] D. Miiller, D. Robaschik, B. Geyer, F.-M. Dittes, J. Horejsi, Fortschr. Phys. 42 (1994) 101. 

[2] X. Ji, Phys. Rev. D 55 (1997) 7114. 

[3] A.V. Radyushkin, Phys. Rev. D 56 (1997) 5524. 

[4] J.C. Collins, L. Frankfurt, M. Strikman, Phys. Rev. D 56 (1997) 2982. 

[5] A.V. Radyushkin, Phys. Lett. B 385 (1996) 333. 

[6] J.C. Collins, A. Freund, Phys. Rev. D 59 (1999) 074009. 

[7] X. Ji, J. Osborne, Phys. Rev. D 58 (1998) 094018. 

[8] F.-M. Dittes, B. Geyer, D. Miiller, D. Robaschik, J. Horejsi, Phys. Lett. B 209 (1988) 325. 

[9] A.V. Efremov, A.V. Radyushkin, Theor. Math. Phys. 42 (1980) 97; Phys. Lett. B 94 (1980) 
245. 

[10] S.J. Brodsky, G.P Lepage, Phys. Lett. B 87 (1979) 359; Phys. Rev. D 22 (1980) 2157. 

[11] M.K. Chase, Nucl. Phys. B 174 (1980) 109. 

[12] Th. Ohrndorf, Nucl. Phys. B 186 (1981) 153. 

[13] V.N. Baier, A.G. Grosin, Nucl. Phys. B 192 (1981) 476. 

[14] A.P Bukhvostov, G.V. Frolov, E.A. Kuraev, L.N. Lipatov, Nucl. Phys. B 258 (1985) 601. 

[15] M.H. Sarmadi, Phys. Lett. B 143 (1984) 471. 

[16] F.-M. Dittes, A.V. Radyushkin, Phys. Lett. B 134 (1984) 359. 

[17] S.V. Mikhailov, A.V. Radyushkin, Nucl. Phys. B 254 (1985) 89. 

[18] L. Frankfurt, A. Freund, V. Guzey, M. Strikman, Phys. Lett. B 418 (1998) 345; (E) ibid. B 
429 (1998) 414. 



[19] I.V. Musatov, A.V. Radyushkin, Evolution and models for skewed parton distributions, |hep- 



ph/9905376 



[20] A.V. Belitsky, B. Geyer, D. Miiller, A. Schafer, Phys. Lett. B 421 (1998) 312. 
[21] L. Mankiewicz, G. Piller, T. Weigl, Eur. J. Phys. C 5 (1998) 119. 

51 



[22] M.K. Chase, Nucl. Phys. B 167 (1980) 125. 

[23] M.A. Shifman, M.I. Vysotsky, Nucl. Phys. B 186 (1981) 475. 

[24] A.G. Shuvaev, Phys. Rev. D 60 (1999) 116005. 

[25] LI. Balitsky, V.M. Braun, Nucl. Phys. B 311 (1989) 541. 

[26] N. Kivel, L. Mankiewicz, Phys. Lett. B 443 (1999) 308; Nucl. Phys. B ??? (1999) ???. 

[27] D. Muller, Phys. Rev. D 58 (1998) 054005. 

[28] A.V. Belitsky, D. Muller, Phys. Lett. B 417 (1997) 129. 

[29] A.V. Belitsky, D. Muller, Nucl. Phys. B 527 (1998) 207. 

[30] D. Muller, Phys. Rev. D 49 (1994) 2525. 

[31] A.V. Belitsky, D. Muller, Nucl. Phys. B 537 (1998) 397. 

[32] A.V. Belitsky, D. Muller, L. Niedermeier, A. Schafer, Phys. Lett. B 437 (1998) 160. 
[33] A.V. Belitsky, D. Muller, L. Niedermeier, A. Schafer, Nucl. Phys. B 546 (1999) 279. 
[34] S.V. Mikhailov, A.V. Radyushkin, Nucl. Phys. B 273 (1986) 297. 
[35] S.V. Mikhailov, Phys. Lett. B 416 (1997) 421. 

[36] B. Geyer, D. Robaschik, M. Bordag, J. Horejsi, Z. Phys. C 26 (1985) 591. 

[37] T. Braunschweig, B. Geyer, J. Horejsi, D. Robaschik, Z. Phys. C 33 (1987) 275. 

[38] S.L. Adler, J.C. Collins, A. Duncan, Phys. Rev. D 15 (1977) 1712. 

[39] J.C. Collins, A. Duncan, S.D. Joglekar, Phys. Rev. D 16 (1977) 438. 

[40] N.K. Nielsen, Nucl. Phys. B 120 (1977) 212. 

[41] K. Symanzik, Comm. Math. Phys. 23 (1971) 49; 
B. Schroer, Nuovo Cim. Lett. 2 (1971) 627; 
G. Parisi, Phys. Lett. B 39 (1972) 643. 

[42] H. Bateman, A. Erdelyi, Higher transcendental functions, V. 1, Mc Graw-Hill, New York 
(1953). 

[43] D. Muller, Z. Phys. C 49 (1991) 293. 

52 



[44] A. V. Belitsky, D. Miiller, Superconformal constraints for QCD conformal anomalies, preprint 
YITP-SB-99-63. 

[45] E.G. Floratos, C. Kounnas, R. Lacaze, Nucl. Phys. B 192 (1981) 417. 

[46] R. Mertig, W.L. van Neerven, Z. Phys. C 70 (1996) 637. 

[47] W. Vogelsang, Phys. Rev. D 54 (1996) 2023; Nucl. Phys. B475 (1996) 47. 

[48] W. Furmanski, R. Petronzio, Phys. Lett. B97 (1980) 437. 

[49] Yu.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

[50] A.V. Belitsky D. Miiller, A. Schafer, Phys. Lett. B 450 (1999) 126. 

[51] G. Curci, W. Furmanski, R. Petronzio, Nucl. Phys. B 175 (1980) 27. 

[52] A.V. Belitsky, A. Freund, D. Miiller, Phys. Lett. B 461 (1999) 270; 
A.V. Belitsky, D. Miiller, Phys. Lett. B 464 (1999) 249. 



53 



